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Abstract 

Motivated  by  recent  efforts  to  mitigate  blast  loading  using  energy-absorbing  materials,  this  paper  uses  analytical  and 
computational  modeling  to  investigate  the  influence  of  mass  distribution  on  the  uniaxial  crushing  of  cellular  sandwich 
plates  under  air  blast  loading.  In  the  analytical  model,  the  cellular  core  is  represented  using  a  rigid,  perfectly-plastic,  lock¬ 
ing  idealization,  as  in  previous  studies,  and  the  front  and  back  faces  are  modeled  as  rigid,  with  pressure  loading  applied  to 
the  front  face  and  the  back-face  unrestrained.  This  model  is  also  applicable  to  the  crushing  of  cellular  media  in  “blast  pen¬ 
dulum”  experiments.  Fluid-structure  interaction  effects  are  treated  using  a  recent  result  that  accounts  for  nonlinear  com¬ 
pressibility  effects  for  intense  air  blasts.  Predictions  of  the  analytical  model  show  excellent  agreement  with  explicit  finite 
element  computations,  and  the  model  is  used  to  investigate  the  response  of  the  system  for  all  possible  distributions  of  mass 
between  the  front  and  back  faces  and  the  cellular  core.  Increasing  the  mass  fraction  in  the  front  face  is  found  to  increase  the 
impulse  required  for  complete  crushing  of  the  cellular  core  but  also  to  produce  undesirable  increases  in  back-face  accel¬ 
erations.  Optimal  mass  distributions  for  mitigating  shock  transmission  through  the  sandwich  plate  are  investigated  by 
maximizing  the  impulse  capacity  while  limiting  the  back-face  accelerations  to  a  specified  level. 

Published  by  Elsevier  Ltd. 
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1.  Introduction 

Sandwich  plates  with  cellular  metal  cores  are  being  widely  considered  for  blast-mitigation  applications,  due 
largely  to  the  energy  absorption  capacity  of  the  cellular  core  material.  Computational  simulations  have  shown 
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that  under  blast  loading,  sandwich  plates  exhibit  reduced  deflections  relative  to  solid  plates  with  the  same  total 
mass  (Xue  and  Hutchinson,  2003,  2004;  Fleck  and  Deshpande,  2004;  Qiu  et  al.,  2004;  Hutchinson  and  Xue, 
2005;  Vaziri  and  Hutchinson,  2006).  These  studies  considered  identical  face  sheets,  and  while  the  influence  of 
varying  the  mass  ratio  between  the  face  sheets  and  the  core  has  been  studied,  the  influence  of  varying  the  mass 
ratio  between  the  two  face  sheets  has  not.  The  influence  of  varying  the  thickness  of  the  face  sheet  nearest  the 
blast,  in  particular,  is  not  obvious,  because  two  competing  effects  are  at  work.  On  the  one  hand,  reducing  the 
thickness  of  this  front  face  means  that  the  blast  impulse  is  imparted  to  a  smaller  mass,  resulting  in  increased 
kinetic  energy,  which  must  be  dissipated  through  crushing  of  the  core.  This  effect  was  noted  by  Xue  and 
Hutchinson  (2004)  in  a  computational  study  of  blast  loading  on  sandwich  plates,  in  which  the  kinetic  energy 
imparted  to  a  sandwich  plate  was  observed  to  be  greater  than  for  a  monolithic  solid  plate  of  the  same  total 
mass.  On  the  other  hand,  reducing  the  front  face  thickness  enhances  the  beneficial  effects  of  fluid-structure 
interaction  (FSI),  which  leads  to  reductions  in  the  impulse  imparted  to  the  system.  A  recent  study  that 
accounted  for  nonlinear  compressibility  effects  in  air  blast  loading  on  freestanding  solid  plates  found  that  such 
impulse  reductions  due  to  FSI  could  be  fairly  significant  (Kambouchev  et  al.,  2006),  more  significant  than  pre¬ 
dicted  by  the  linear  acoustic  FSI  analysis  of  Taylor  (1963),  which  has  been  used  in  previous  analyses  of  blast 
loading  on  sandwich  plates.  Depending  which  of  these  effects  is  dominant,  benefits  might  be  achieved  by  either 
increasing  or  decreasing  the  mass  fraction  in  the  face  sheet  nearest  the  blast. 

Motivated  by  these  observations,  this  paper  uses  analytical  and  computational  modeling  to  investigate  the 
influence  of  mass  distribution  on  the  uniaxial  crushing  of  sandwich  plates  under  air  blast,  considering  all  pos¬ 
sible  distributions  of  mass  between  the  two  face  sheets  and  the  cellular  core.  While  previous  studies  of  blast 
loading  on  sandwich  plates  have  primarily  focused  on  minimizing  dynamic  deflections,  the  main  interest  in  the 
present  study  is  in  mitigating  shock  transmission  through  the  thickness  of  the  sandwich  plate.  A  potential 
application  is  for  shock  isolation  of  electronic  components  mounted  to  the  back  face  of  an  armor  panel.  Accel¬ 
erations  of  the  back  face  are  of  particular  concern  in  such  applications,  because  damage  to  an  electronic  com¬ 
ponent  is  related  to  the  acceleration  transmitted  to  its  base. 

Because  stress-strain  curves  for  cellular  media  exhibit  a  characteristic  plateau,  it  has  been  noted  (Hanssen 
et  al.,  2002)  that  a  cellular  layer  could  potentially  serve  a  protective  function  against  blast  loading  by  limiting 
the  stresses  transmitted  through  the  layer  to  the  plateau  stress  erP.  However,  this  protective  function  is  lost  if 
complete  densification  of  the  cellular  layer  is  permitted,  in  which  case  the  transmitted  stresses  can  actually  be 
amplified.  Such  amplification  can  occur  because  the  crushing  mechanism  of  cellular  media  under  blast  or  high¬ 
speed  impact  involves  propagation  of  a  densification  front  with  the  characteristics  of  a  shock  wave.  While  the 
stress  ahead  of  the  densification  front  is  limited  to  the  plateau  stress,  a  velocity-dependent  jump  in  stress  across 
the  densification  front  can  result  much  higher  stress  in  the  crushed  material  behind  the  densification  front. 
This  stress  enhancement  due  to  shock  formation  has  been  exploited  in  the  use  of  metal  foam  projectiles  to 
simulate  shock  loading  on  structures  (Rathbun  et  al.,  2006;  Radford  et  al.,  2006).  Under  blast  loading,  stresses 
several  times  larger  than  the  peak  applied  blast  pressure  can  be  transmitted  through  a  cellular  layer  if  the 
thickness  of  the  layer  is  insufficient  to  arrest  the  densification  front,  as  observed  by  Li  and  Meng  (2002). 

This  mechanism  of  stress  enhancement  due  to  shock  formation  helps  to  account  for  the  surprising  obser¬ 
vation  that  the  use  of  cellular  materials  under  blast  loading  has  in  many  instances  led  to  enhancement,  rather 
than  mitigation,  of  blast  effects.  Cooper  et  al.  (1991)  reported  enhanced  thoracic  injury  in  animals  whose  chest 
cavities  were  “protected”  by  a  foam  layer,  and  computational  studies  conducted  by  Friend  (2005)  revealed 
that  thoracic  injury  is  exacerbated  by  shock  wave  compression  of  alternating  soft  foam  layers  and  subsequent 
acceleration  of  more  rigid  layers  in  porous  foam  body  armor.  Experiments  by  Hanssen  et  al.  (2002)  showed 
that  increased  upswing  resulted  from  the  addition  of  an  aluminum  foam  layer  to  the  face  of  a  massive  “pen¬ 
dulum”  subjected  to  blast  loading,  and  other  examples  of  enhanced  blast  effects  are  described  by  Li  and  Meng 
(2002).  In  discussing  the  blast  pendulum  experiments  of  Hanssen  et  al.  (2002),  Nesterenko  (2002)  credited  the 
observed  impulse  enhancement  to  the  first  of  the  two  competing  effects  mentioned  above,  noting  that  the  blast 
impulse  would  have  been  imparted  primarily  to  a  lightweight  plate  covering  the  foam  layer,  leading  to  signif¬ 
icantly  higher  kinetic  energy  than  if  the  same  impulse  were  imparted  directly  to  the  more  massive  pendulum. 
Nesterenko  (2002)  suggested  that  complete  crushing  could  then  lead  to  detachment  and  rebounding  of  the 
densified  foam  layer,  thus  increasing  the  resultant  momentum  of  the  pendulum.  It  is  noted  that  the  blast  pen¬ 
dulum  can  be  considered  an  extreme  example  of  a  “sandwich  plate”  with  unequal  face  sheets,  with  most  of  the 
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mass  in  the  back  face  (i.e.,  the  pendulum)  and  relatively  small  masses  in  the  front  face  and  the  cellular  core. 
The  present  study  can  thus  afford  insights  into  the  influence  of  mass  distribution  on  the  poor  performance  of 
this  system  under  blast  loading. 

To  investigate  the  influence  of  mass  distribution  on  the  uniaxial  crushing  of  sandwich  plates  under  air  blast 
loading,  an  analytical  model  is  developed  in  this  paper  in  which  the  front  and  back  faces  are  treated  as  rigid, 
with  a  pressure  pulse  applied  to  the  front  face  and  the  back-face  unrestrained,  and  the  cellular  core  is  repre¬ 
sented  using  the  well-established  rigid,  perfectly-plastic,  locking  (R-P-P-L)  continuum  idealization.  The  R-P- 
P-L  idealization  was  originally  developed  to  model  the  dynamic  crushing  of  wood  (Reid  and  Peng,  1997)  and 
has  subsequently  been  applied  to  metal  honeycombs  and  foams  in  a  number  of  studies  (Harrigan  et  ah,  1999; 
Hanssen  et  ah,  2002;  Tan  et  ah,  2005b).  The  R-P-P-L  idealization  can  be  applied  to  a  variety  of  cellular  topol¬ 
ogies,  provided  that  the  cell  dimensions  are  sufficiently  small  relative  to  the  domain  of  interest  to  justify  the 
continuum  approximation.  Based  on  impact  experiments  for  aluminum  foams,  Tan  et  ah  (2005a)  observed 
that  the  domain  of  interest  should  include  at  least  five  average  cell  diameters  for  a  valid  continuum  approx¬ 
imation.  For  cases  in  which  the  cell  sizes  are  large  and  the  individual  core  elements  are  relatively  stocky,  rep¬ 
resenting  the  core  crushing  as  a  planar  densification  front  propagating  through  a  continuum  may  be 
inappropriate,  and  the  core  crushing  may  be  better  represented  using  the  model  developed  by  Vaughn 
et  ah  (2005),  which  considers  coupled  plastic  wave  propagation  and  nonlinear  buckling  of  individual  columns. 
Xue  and  Hutchinson  (2006)  showed  that  the  model  of  Vaughn  et  ah  (2005)  predicted  quite  accurately  the 
dynamic  crushing  strength  of  square  honeycomb  sandwich  cores  under  different  crushing  rates.  While  the 
R-P-P-L  idealization  does  not  explicitly  incorporate  rate  dependence  in  the  plateau  stress,  rate-dependent 
behavior  results  from  the  shock  wave  formation  and  the  associated  velocity-dependent  stress  jump  mentioned 
previously.  Tan  et  ah  (2005a,  Fig.  18),  showed  that  a  one-dimensional  shock  model  using  the  R-P-P-L  ideal¬ 
ization  was  able  to  accurately  predict  experimental  observations  of  increasing  crushing  stress  with  increasing 
impact  velocity  for  aluminum  foam  projectiles. 

As  illustrated  in  Fig.  1,  a  number  of  previous  studies  have  developed  analytical  models  to  represent  the  one¬ 
dimensional  crushing  of  cellular  media  under  various  dynamic  loading  conditions.  Reid  and  Peng  (1997)  used 
the  R-P-P-L  idealization  to  model  the  dynamic  crushing  of  wood  cylinders  impacting  a  rigid  wall,  obtaining  an 
expression  for  the  enhancement  of  stress  transmitted  to  the  rigid  wall  as  a  consequence  of  the  stress  jump 
across  the  shock  front.  This  model  was  subsequently  applied  to  aluminum  honeycomb  projectiles  by  Harrigan 
et  ah  (1999).  Hanssen  et  ah  (2002)  investigated  the  crushing  of  a  metal  foam  bar  with  a  triangular  pressure 
pulse  applied  to  a  rigid  cover  sheet,  comparing  finite  element  results  with  an  analytical  solution  using  the 
R-P-P-L  idealization.  They  obtained  an  expression  for  the  minimum  thickness  required  to  fully  absorb  a  given 
impulse,  thereby  limiting  the  stress  transmitted  to  the  rigid  support  to  the  plateau  stress  erP.  Li  and  Meng 
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Fig.  1.  One-dimensional  systems  considered  in  previous  studies. 
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(2002)  used  a  nonlinear  spring-mass  model  to  examine  the  shock  wave  propagation  through  a  cellular  medium 
subject  to  a  rectangular  pressure  pulse,  and  as  noted  previously,  they  observed  a  sharp  increase  in  the  stress 
transmitted  to  the  rigid  support  for  impulses  exceeding  the  energy  absorption  capacity  of  the  cellular  layer. 
The  R-P-P-L  idealization  was  extended  by  Lopatnikov  et  al.,  2003,  2004,  who  used  an  elastic,  perfectly-plastic, 
rigid  (E-P-P-R)  idealization  to  model  the  dynamic  deformation  of  metal  foam  materials  under  both  Taylor 
impact  and  plate  impact  boundary  conditions.  Li  and  Reid  (2006)  subsequently  proposed  corrections  to 
the  unloading  assumption  used  in  the  E-P-P-R  model  and  the  impact  velocity  regimes  defined  by  Lopatnikov 
et  al.  (2003)  for  the  plate  impact  problem. 

In  analyzing  the  dynamic  response  of  sandwich  beams  and  plates  to  blast  loading,  Fleck  and  Deshpande 
(2004)  and  Hutchinson  and  Xue  (2005)  divided  their  analyses  into  three  temporal  stages,  of  which  stage  I 
involved  FSI  in  the  reflection  of  the  incident  pressure  pulse  from  the  front  face  of  the  sandwich  panel,  and 
stage  II  involved  one-dimensional  compaction  of  the  cellular  core  with  initial  velocities  resulting  from  the 
stage  I  analysis.  Both  studies  treated  FSI  in  stage  I  by  applying  the  linear  acoustic  analysis  of  Taylor 
(1963)  to  the  front  face,  while  Hutchinson  and  Xue  (2005)  extended  this  analysis  to  incorporate  the  resistance 
provided  by  a  perfectly-plastic  core.  Both  studies  analyzed  the  core  crushing  in  stage  II  by  employing  a  quasi¬ 
static  assumption  for  the  energy  dissipation  in  the  core,  with  Fleck  and  Deshpande  (2004)  employing  the  R-P- 
P-L  idealization  and  Hutchinson  and  Xue  employing  a  perfectly-plastic  idealization  with  no  compaction  limit. 
In  a  subsequent  analysis  of  shock  wave  propagation  in  cellular  media,  however,  Radford  et  al.  (2005)  noted 
that  the  energy  dissipation  in  the  core  could  be  significantly  larger  than  the  quasi-static  value  because  of  addi¬ 
tional  dissipation  within  the  shock  wave.  This  additional  dissipation  is  associated  with  a  drop  in  kinetic  energy 
across  the  shock  front,  which  follows  from  the  balance  of  momentum,  not  from  an  explicit  assumption  of 
material  rate  dependence.  Assuming  a  linear  rate  dependence  for  the  cellular  material,  Radford  et  al. 
(2005)  presented  an  asymptotic  analysis  of  the  shock  structure,  which  enabled  a  quantification  of  the  viscous 
dissipation  that  balances  the  drop  in  kinetic  energy.  Tan  et  al.  (2005b)  also  noted  the  incorrectness  of  a  global 
energy  balance  that  assumes  quasi-static  energy  dissipation  in  the  core,  using  a  thermo-mechanical  approach 
to  derive  the  governing  equations  for  shock  propagation  in  an  R-P-P-L  solid  and  showing  that  quasi-static 
assumption  leads  to  a  contradiction.  The  velocity-dependent  energy  dissipation  associated  with  the  shock 
front  was  properly  accounted  for  in  a  more  recent  model  developed  by  Deshpande  and  Fleck  (2005)  for 
one-dimensional  crushing  of  sandwich  plates.  This  model  also  relaxed  the  requirement  of  a  strict  temporal  sep¬ 
aration  between  FSI  in  stage  I  and  core  crushing  in  stage  II  by  considering  time -varying  loading  on  the  front 
face  rather  than  simply  applying  the  resulting  impulse  to  the  front  face  as  an  initial  velocity. 

In  terms  of  the  representation  of  shock  wave  propagation  and  associated  energy  dissipation  in  the  cellular 
core,  the  present  model  is  equivalent  to  that  of  Deshpande  and  Fleck  (2005),  although  it  was  developed  inde¬ 
pendently  and  was  previously  reported  briefly  in  Main  and  Gazonas  (2005).  Key  differences  are  that  Deshpan¬ 
de  and  Fleck  (2005)  considered  face  sheets  of  equal  mass,  while  the  present  model  allows  for  arbitrary  masses 
of  the  front  and  back  faces,  and  that  Deshpande  and  Fleck  (2005)  used  the  linear  acoustic  analysis  of  Taylor 
(1963)  to  represent  FSI  effects  for  underwater  shock  loading,  while  the  present  model  incorporates  the  recent 
analysis  of  Kambouchev  et  al.  (2006)  to  account  for  nonlinear  compressibility  and  FSI  effects  in  air-blast  load¬ 
ing.  Vaziri  and  Hutchinson  (2006)  presented  two  approaches  for  implementing  the  result  of  Kambouchev  et  al. 
(2006)  in  computational  simulations  of  sandwich  plates  under  air-blast  loading,  one  in  which  the  reflected 
impulse  is  imparted  to  the  front  face  through  an  initial  velocity  and  another  in  which  a  time-dependent  pres¬ 
sure  history  is  applied  to  the  front  face.  The  present  study  uses  a  modification  of  the  latter  approach  that 
accounts  for  the  initial  resistance  provided  by  the  plateau  stress  of  the  cellular  core. 

The  outline  of  this  paper  is  as  follows.  The  derivation  of  the  analytical  model  is  first  presented  in  Section  2, 
and  in  Section  3  predictions  of  the  analytical  model  are  compared  with  finite  element  results.  Initiation  and 
arrest  of  densification  are  considered  in  Section  4,  in  which  the  analytical  model  is  used  to  investigate  the  influ¬ 
ence  of  mass  distribution  on  the  critical  reflected  impulse  required  for  complete  densification  of  the  cellular 
core  under  different  levels  of  peak  reflected  pressure.  The  influence  of  mass  distribution  on  the  time  required 
for  complete  densification  and  the  corresponding  energy  dissipation  in  the  core  are  also  investigated.  FSI 
effects  are  considered  in  Section  5,  in  which  the  influence  of  mass  distribution  on  the  critical  incident  impulse 
is  investigated  and  the  work  done  by  the  air  blast  on  the  sandwich  plate  is  compared  with  the  work  done  on  a 
corresponding  solid  plate  with  the  same  total  mass.  Finally,  accelerations  of  the  back  face  are  investigated  in 
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Section  6,  and  optimal  mass  distributions  for  mitigation  of  shock  transmission  are  investigated  by  seeking  to 
maximize  the  impulse  that  can  be  absorbed  while  limiting  the  back-face  accelerations  to  a  specified  level. 

2.  Analytical  model 

A  strip  through  the  thickness  of  a  sandwich  plate  is  considered,  as  depicted  with  a  partially  compacted  core 
in  Fig.  2a.  The  R-P-P-L  idealization,  depicted  in  Fig.  2b,  is  used  to  represent  the  stress-strain  relationship  of 
the  core  material.  The  applied  pressure  pulse  pr(  t)  is  the  pressure  reflected  from  the  front  face  of  the  sandwich 
panel  and  is  influenced  by  the  nonlinear  compressibility  of  air  and  by  fluid-structure  interaction  (FSI). 
Because  of  FSI  effects,  pr  actually  depends  on  the  motion  of  the  front  face  and  cannot  be  formally  defined 
as  a  function  solely  of  time  t.  However,  Section  5  presents  an  approximate  approach,  based  on  the  analysis 
of  Kambouchev  et  al.  (2006),  whereby  nonlinear  compressibility  and  FSI  effects  are  accounted  for  in  defining 
a  reflected  pressure  pulse  pT(t)  that  results  from  a  specified  incident  pressure  pulse  pf).  While  no  particular 
functional  form  is  assumed  for  pr{  t)  in  the  derivation  of  the  analytical  model,  it  is  assumed  that  the  peak 
reflected  pressure  pR  occurs  initially,  so  that  /MO)  =  pR  and  pT(t)  M  pK  for  t  >  0. 

2.1.  Equation  of  motion 

A  strip  of  sandwich  plate  with  unit  cross-sectional  area  is  considered  (i.e.,  with  unit  area  exposed  to  the 
blast  pressure),  with  total  mass  given  by  m  =  mx  +  p0£0  +  m2,  where  p0  and  £0  are  the  uncompressed  density 
and  initial  thickness  of  the  cellular  core,  and  m  i  and  m2  are  the  areal  densities  of  the  front  and  back  faces, 
respectively.  The  acceleration  of  the  center  of  mass,  denoted  wG,  follows  directly  from  application  of  Newton’s 
second  law  to  the  entire  system,  whereby  all  internal  forces  cancel  and  the  time  derivative  of  the  net  momen¬ 
tum  is  equivalent  to  the  net  external  force  (because  a  strip  of  unit  cross-sectional  area  is  considered,  forces  and 
pressures  can  be  used  interchangeably): 

pft)=mua  (1) 

The  total  momentum  of  the  system  can  be  obtained  by  integrating  Eq.  (1)  with  respect  to  time,  assuming  that 
the  panel  is  initially  at  rest,  wG( 0)  =  0,  to  give  the  impulse-momentum  relationship 

zr(f)  =  mitQ  (2) 

where  /,( t)  denotes  the  reflected  impulse  per  unit  area: 

h W  =  [  Pr{i)di  (3) 

Provided  that  the  applied  pressure  is  sufficiently  high,  densification  of  the  cellular  core  commences  at  the 
front  face,  and  a  shock  front  propagates  through  the  core.  (The  pressure  required  to  initiate  densification 
is  investigated  subsequently.)  According  to  the  simplified  model  of  Fig.  2b,  the  densified  core  material  moves 
as  a  rigid  body  with  the  same  velocity  as  the  front  face,  while  the  uncompressed  core  material  moves  as  a  rigid 
body  with  the  velocity  of  the  back  face.  Since  the  motion  of  the  center  of  mass  is  known  from  Eqs.  (1)  and  (2), 
it  is  convenient  to  express  the  positions  of  the  front  and  back  faces  of  the  panel  with  reference  to  the  position 
of  the  mass  center,  as  follows: 


Fig.  2.  Analytical  model  definition:  (a)  strip  of  sandwich  panel  with  partially  compacted  core  and  (b)  engineering  stress-strain  relationship 
for  cellular  core  material  (R-P-P-L  idealization). 
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mi  =  mg  —  (£  —  xg);  »2  =  uq  +xq  (4a,  b) 

where  m,  and  m2  denote  the  positions  of  the  front  and  back  faces,  respectively,  £  denotes  the  total  thickness  of 
the  core,  and  xG  denotes  the  distance  from  the  center  of  mass  to  the  back  face,  as  shown  in  Fig.  2a.  With  sD 
denoting  the  volumetric  strain  at  complete  densification,  it  follows  from  conservation  of  mass  that  the  density 
of  the  compressed  core  material  is  p0/{  1  —  £d)  and  the  total  thickness  of  the  core  is  given  by 

£  =  £0  -  sD(£o  -  x)  (5) 

where  x  denotes  the  thickness  of  the  uncompressed  core  material,  and  the  thickness  of  the  shock  front  itself  is 
assumed  to  be  negligible.  By  forming  and  differentiating  an  expression  for  xG  the  following  relation  can  be 
obtained: 

xg  =  (,£T>/m)[r>h  +  Po{£o  ~  x)]x  (6) 

The  following  expressions  for  the  velocities  of  the  front  and  back  faces  can  then  be  obtained  by  differentiating 
Eqs.  (4a,b)  and  (5)  and  combining  with  Eq.  (6): 

mi  =  [ir(t)  -  £d(Pox  +  m2)x]/m\  m2  =  {4(0  +  £d[«i  +  Po(£o  ~  x)]x}/m  (7a,  b) 

in  which  the  substitution  mg  =  iT{t)/m  from  Eq.  (2)  has  been  made. 

The  stress  just  ahead  of  the  shock  front  is  erP,  and  application  of  Newton’s  second  law  to  the  material  ahead 
of  the  shock  front  yields  the  following  equation: 

op  =  (p0x  +  m2)u2  (8) 

Differentiating  Eq.  (7b)  and  substituting  into  Eq.  (8)  then  yields  the  following  nonlinear  ordinary  differential 
equation  for  x: 

— £D[»zi  +  p0(£o  -  x)]x  +  £Dp0x2  =  pr(t)  -  oPm/ (p0x  +  m2)  (9) 

Eq.  (9)  can  be  integrated  numerically  with  initial  conditions  x(0)  =  £0  and  x(0)  =  0. 


2.2.  Energy  balance 

The  balance  of  mechanical  energy  for  the  system  of  Fig.  2  requires  that  the  work  done  by  the  pressure  pulse, 
denoted  Wp,  must  equal  the  sum  of  the  kinetic  energy  of  the  system,  denoted  Ek,  and  the  energy  dissipated  in 
the  core  material,  denoted  Ed: 


Wp  =  Ek  +  Ed  (10) 

The  work  per  unit  area  done  by  the  pressure  pulse  can  be  expressed  as  follows,  in  terms  of  the  applied  pres¬ 
sure  pT(t)  and  the  velocity  of  the  front  face  up. 

Wp  =  f  pT(i)ui  dt  (11) 

Jo 

The  expression  for  m,  from  Eq.  (7b)  can  then  be  substituted  into  Eq.  (11).  The  resulting  integral  contains  the 
term  iT(t)pr(t),  which  integrates  directly  to  yielding  the  following  expression: 


1  r< 

WP  =  ^‘riO/m  -  (£d/m)  J  pr(t)(PoX  +  m2)j 


d  t 


(12) 


The  first  term  in  Eq.  (12)  represents  the  work  done  by  the  pressure  pulse  if  no  densification  occurs  (x  =  0),  in 
which  case  the  front  face  moves  at  the  same  velocity  as  the  center  of  mass,  u\  =  uQ.  The  second  term  in  Eq. 
( 12)  represents  the  enhancement  of  this  work  resulting  from  densification  of  the  cellular  core.  This  term  is  po¬ 
sitive  because  x  <  0  as  the  shock  front  propagates.  Densification  of  the  core  results  in  larger  velocities  of  the 
front  face  (ii\  >  mg),  so  that  the  applied  pressure  acts  through  a  larger  distance,  resulting  in  increased  work. 
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With  the  thickness  of  the  shock  front  assumed  negligible,  the  kinetic  energy  of  the  sandwich  panel  per  unit 
area  can  be  expressed  as  follows,  as  the  sum  of  contributions  from  the  portions  on  each  side  of  the  shock  front: 

1  1 

£k  =  ^  imi  +  Po^o  -  x)]u-{  + -(p0x  +  m2)u2  (13) 

Substituting  the  expressions  for  ii\  and  u2  from  Eqs.  (7a, b)  and  collecting  terms  yields  the  following  expression 
for  the  kinetic  energy: 


Ek  =  1 'r  (t)/m  +  *  (el/m)[mi  +  p0(£0  -  x)\(p0x  +  m2)i 


(14) 


The  first  term  in  Eq.  ( 14)  represents  the  kinetic  energy  associated  with  motion  of  the  mass  center  and  is  equiv¬ 
alent  to  the  first  term  in  Eq.  (12).  The  second  term  in  Eq.  (14)  is  zero  if  no  densification  occurs  (x  =  0),  and  it 
represents  an  enhancement  of  the  kinetic  energy  resulting  from  relative  motion  of  the  front  and  back  faces. 

Because  the  uncompressed  material  ahead  of  the  shock  front  and  the  densified  material  behind  the  shock 
front  are  both  assumed  to  act  rigidly,  it  is  only  within  the  shock  front  that  energy  dissipation  occurs.  An 
expression  for  the  energy  dissipated  in  the  core  can  then  be  obtained  by  considering  a  control  volume  that 
encloses  the  shock  front.  With  the  thickness  of  the  shock  front  assumed  negligible,  no  mass  or  kinetic  energy 
can  accumulate  within  the  control  volume.  Therefore,  the  flux  of  kinetic  energy  into  the  control  volume  must 
be  balanced  by  the  rate  of  external  work  on  the  control  volume  and  the  rate  of  energy  dissipation  within  the 
control  volume.  The  rate  of  energy  dissipation  within  the  control  volume  can  then  be  expressed  as  follows: 

Ed  =  -  iii) (iq  ~  «i)/£d  -  +  (ffp  +  Acr)w!  (15) 


The  first  term  on  the  right-hand  side  gives  the  flux  of  kinetic  energy  into  the  control  volume,  where 
p0(ui  -  u2)/et>  gives  the  flux  of  mass  through  the  control  volume,  \u\  represents  the  kinetic  energy  per  unit 
mass  entering  the  control  volume,  and  (  u\  represents  the  kinetic  energy  per  unit  mass  leaving  the  control  vol¬ 
ume.  The  remaining  terms  give  the  net  rate  of  external  work  on  the  system  within  the  control  volume.  The 
quantity  Act  represents  the  jump  in  stress  across  the  shock  front,  and  it  follows  from  conservation  of  momen¬ 
tum  that  An  is  given  by 

Act  =  P0(m2  —  «i)2/ed  (16) 


Substituting  Eq.  (16)  into  Eq.  (15)  and  collecting  terms  yields  the  following  expression: 

1  3 

E d  =  -Op(w2  -  til)  -  ^Po(ui  -  Ml)  /eD 


(17) 


Noting  that  u2  —  it\  =  t,  it  follows  from  differentiation  of  Eq.  (5)  that  u2  —  ii\  =  eDx,  whereby  Eq.  (17)  can  be 
expressed  in  terms  of  x  as 

1 

Ed  =  -GpSox  -  ^Poedx  (18) 


With  the  conditions  Ed  =  0  and  v;  =  £0  at  t  =  0,  Eq.  ( 18)  can  then  be  integrated  to  give 


Ed  —  o>£d(^o  ~  x)  —  -  PoeD 


i3  d  t 


(19) 


The  first  term  in  Eq.  (19)  represents  the  quasi-static  energy  dissipation  associated  with  densification  of  the  core 
under  a  constant  stress  erP,  while  the  second  term  represents  a  rate-dependent  enhancement  of  energy  dissipa¬ 
tion  associated  with  the  shock  front. 

It  is  noted  that  the  equation  of  motion  (9)  can  be  derived  in  an  alternative  manner  by  using  the  work  and 
energy  expressions  presented  in  this  section.  Differentiation  of  Eq.  (10)  with  respect  to  time  yields  the  follow¬ 
ing  relationship: 

W  p  =  E k  +  Ed 


(20) 
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It  can  then  be  readily  verified  that  differentiation  of  Eqs.  (12),  (14),  and  (19)  and  substitution  into  Eq.  (20) 
yields  the  same  differential  equation  as  Eq.  (9).  The  equation  of  motion  (9)  is  therefore  fully  consistent  with 
the  mechanical  energy  balance  presented  in  this  section. 

2. 3.  Nondimensionalization 


For  the  sake  of  generality,  it  is  useful  to  express  the  equation  of  motion  (9)  in  terms  of  nondimensional 
quantities.  A  nondimensional  coordinate  £  is  introduced  to  denote  the  position  of  the  shock  front: 

Z  =  x/eo  (2i) 

It  is  noted  that  £  represents  the  remaining  mass  fraction  of  uncompressed  core  material  and  that  £  =  1  ini¬ 
tially.  The  symbol  zR  is  introduced  to  denote  the  total  reflected  impulse  per  unit  area: 


so  that  4(0  — >  /r  as  I  — >  oo.  A  nondimensional  time  t  can  be  defined  as  follows: 
t  =  to? /iR 


(22) 

(23) 


The  derivatives  of  the  shock  front  position  x  can  then  be  expressed  in  terms  of  nondimensional  variables  as 
x  =  (ioOp/ip/jd^/dx  and  x  =  (£0ol  /  i^)d2  ^  /  dx2 .  A  nondimensional  reflected  pressure  Pr(x)  can  also  be  defined 
as  a  function  of  nondimensional  time  through  division  of  /;,-( t)  by  the  plateau  stress  erP  and  a  change  of  vari¬ 
ables  using  Eq.  (23): 


^r(t) 


pA  0 


(Tp 


f=TiR/uP 


(24) 


For  convenience  in  numerical  solution,  the  additional  state  variable  v  —  d^/dr  is  introduced,  which  represents 
the  nondimensional  relative  velocity  of  the  shock  front,  whereby  the  equation  of  motion  can  be  recast  as  a 
first-order  differential  equation  in  vector  form  as  follows: 


d^ 

dt 


£ 

v 


-hjPrbH’ioi+'k)  ‘l+'Jp"2 

loii+ijjh-f) 


(25) 


where  /R  is  a  nondimensional  total  impulse,  defined  as 


!r  /  Po 
m  Y  o>£d 


(26) 


and  rj0,  i]\.  and  t]2  represent  the  fractions  of  the  total  mass  in  the  core,  the  front  face,  and  the  back  face, 
respectively: 

>lo  =  Po^o/m',  rh=ml/m\  rj2=m2/m  (27) 

For  a  given  form  of  the  pressure  input  Pr{x),  Eq.  (25)  can  be  integrated  numerically  with  the  initial  conditions 
£(0)  =  1  and  u(0)  =  0,  corresponding  to  initially  uncrushed  foam. 

Once  time  histories  of  £  and  v  have  been  obtained  by  numerical  integration,  it  is  necessary  to  express  the 
response  quantities  of  interest  in  terms  of  these  nondimensional  state  variables  and  the  nondimensional  pres¬ 
sure  Pr( r).  Nondimensional  front-face  and  back-face  velocities  can  be  defined  as  follows: 

hi  h2 

Vi  =  /  ;  v2  = 

i  r/'«  ir /m 

where  ipjm  is  the  limiting  velocity  of  the  center  of  mass  as  t  — >  oo,  which  follows  from  Eqs 
substituting  Eqs.  (7a, b)  into  Eqs.  (28a, b)  the  following  expressions  can  be  obtained: 

g,=  ;,(,)  +  ''»(■  -«  +  W.„ 


(28a,  b) 
.  (2)  and  (22).  By 


h  =4W 


(29a,  b) 
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where  /,(t)  is  a  nondimensional  impulse  defined  as  follows,  so  that  /r(t)  — >  1  as  t  — >  oo: 


/r(t)  =  [  PT( i)dt 
Jo 


h(0 


t=TiK/oP 


(30) 


The  acceleration  of  the  back-face  ii2  is  of  interest  in  evaluating  the  protective  function  of  the  cellular  core, 
and  it  follows  from  Eq.  (8)  that  this  is  given  by  ii2  =  oo/(p0x  +  m2).  A  nondimensional  back-face  acceleration 
can  then  be  defined  as  follows: 


a  2  = 


«2 


1 


cjp/m  t]  0£  +  rj  2 


(31) 


Nondimensional  versions  of  the  work  and  energy  quantities  from  Section  2.2  can  be  defined  through  divi¬ 


sion  by  i\/(lm),  which  represents  the  kinetic  energy  associated  with  the  center  of  mass  in  the  limit  as  t 


oo: 


W  n  = 


Wn 


2  R 


A 


£k  = 


ArA 


ea  = 


1  ;2 


?rA 


(32a,  b,  c) 


Substituting  Eq.  ( 14)  into  Eq.  (32b),  the  nondimensional  kinetic  energy  can  be  expressed  as  follows: 

Ek  =  AW  +  Ao/4)['h  +  'M1  -  C)](»/2  +  '7oA2  (33) 


Because  the  expressions  for  Wp  and  Ed  in  Eqs.  (12)  and  (19)  contain  integrals  with  respect  to  time,  it  is  con¬ 
venient  to  express  the  nondimensional  quantities  Wp  and  Ed  in  terms  of  their  derivatives  with  respect  to  non- 
dimensional  time  t,  as  follows: 

d  W 

=  2P,(t)EMt)  -  K/4)(^  +  r,2)v]  (34) 

^  =  -(2)/0/4)u  -  ('i2/4)°3  (35) 

In  numerical  solution  for  Wp  and  Ed,  the  state  vector  in  Eq.  (25)  can  then  be  augmented  to  include  Wp  and  Ed, 
whereby  Eqs.  (25),  (34),  and  (35)  can  be  integrated  simultaneously.  The  following  expression  for  the  nondi¬ 
mensional  quasi-static  energy  dissipation,  corresponding  to  the  first  term  in  Eq.  (19),  can  be  obtained  by  inte¬ 
grating  the  first  term  in  Eq.  (35): 

ET  =  (2%/4)(l  -  9  (36) 

A  nondimensional  peak  reflected  pressure  PR  can  be  defined  as 

Pr=Pr/gp  (37) 


where  pR  =  /;,(())  is  the  peak  reflected  pressure,  assumed  to  occur  initially.  While  the  derivations  thus  far  have 
not  assumed  any  particular  functional  form  for  the  reflected  pressure  pulse  pT(  t),  the  computational  simula¬ 
tions  presented  in  Section  3  consider  a  triangular  pressure  pulse  decreasing  linearly  from  the  peak  pressure 
pR  at  t  —  0  to  zero  at  t  =  tR.  Using  Eqs.  (24)  and  (30),  the  nondimensional  pressure  and  impulse  for  the  tri¬ 
angular  pulse  can  be  expressed  as  follows: 


J  PR(1  -  t/tr),  0  <  t  ^  tr  r  (t/tr)(2  -  t/tR),  0^t<tr 

\  o,  T  >  Tr  l  1,  r  Si  Tr 


(38a,  b) 


where  the  nondimensional  impulse  duration  is  given  by  rR  =  ?RerP//R,  according  to  Eq.  (23),  and  because 
zR  =pRrR/ 2  for  the  triangular  pulse,  it  follows  that  rR  =  2 /PR. 

The  approximate  treatment  of  FSI  in  Section  5  considers  an  exponential  reflected  pressure  pulse  of  the  form 
pr(t)  =  pRe~‘/lR,  where  ?R  is  the  decay  period.  Using  Eqs.  (24)  and  (30),  the  nondimensional  pressure  and 
impulse  for  the  exponential  pulse  can  be  expressed  as  follows: 

Pr( t)  =  PR  e~T/TR;  7r(t)  =  1  -  e-z/XR 

where  tr  =  1  /PR  according  to  Eq.  (23)  with  zR  =  pRtR  for  the  exponential  pulse. 


(39a,  b) 
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The  limiting  case  of  a  zero  duration  impulse  can  be  considered  by  setting  pv(t)  =  iRd{  t)  and  4(0  =  Ir H(t), 
where  S(t)  is  the  Dirac  delta  function  and  H(t)  is  the  Heaviside  step  function.  The  nondimensional  pressure 
and  impulse  in  this  limit  are  simply  given  by  Pr(x)  =  <5(t)  and  IT(x)  =  //( t),  respectively.  By  integrating  the 
equation  of  motion  (25)  about  x  =  0,  it  can  be  shown  that  the  Dirac  delta  pressure  pulse  imparts  a  nondi¬ 
mensional  relative  velocity  of  v0  =  —A/^h'la)  to  the  shock  front.  The  equation  of  motion  (25)  in  this  lim¬ 
iting  case  can  then  be  solved  by  imposing  the  initial  condition  n(0)  =  v0  and  setting  P,-(r)  =  0. 

3.  Comparison  with  computational  predictions 

In  this  section,  the  predictions  of  the  analytical  model  are  compared  with  explicit  finite  element  com¬ 
putations  using  LS-DYNA.  The  computational  simulations  follow  fairly  closely  those  presented  by  Hans- 
sen  et  al.  (2002),  although  the  present  study  considers  an  unrestrained  back  face  with  different  values  of 
the  front-face  and  back-face  masses,  while  Hanssen  et  al.  considered  a  fixed  back  face.  In  the  computa¬ 
tions,  the  cellular  core  was  represented  by  a  single  row  of  solid  elements  with  total  thickness  £o  =  5  cm, 
using  material  model  26  (*MAT_HONEYCOMB)  with  p0  =  250  kg/m3,  <xP  =  1  MPa,  and  e0  =  0.7.  A  large 
elastic  modulus  of  E  =  700  GPa  was  used  to  represent  the  “rigid”  portions  of  the  idealized  stress-strain 
relationship  in  Fig.  2b,  and  Poisson’s  ratio  was  set  to  zero.  The  solid  elements  were  defined  with  an  initial 
aspect  ratio  of  2:1  (axial  to  transverse),  resulting  in  an  aspect  ratio  of  0.6:1  at  complete  densification.  The 
front  and  back  faces  were  represented  in  the  computations  by  added  nodal  masses,  and  two  different  mass 
distributions  were  considered,  as  indicated  in  Table  1.  The  “blast  pendulum”  case  corresponds  to  the  blast 
pendulum  experiments  of  Hanssen  et  al.  (2002),  with  the  large  back-face  mass  representing  the  pendulum. 
The  “sandwich  plate”  case,  with  equal  front-face  and  back-face  masses,  corresponds  to  a  sandwich  plate 
as  considered  by  Xue  and  Hutchinson  (2003,  2004)  and  Fleck  and  Deshpande  (2004). 

FSI  effects  were  not  explicitly  considered  in  the  computational  simulations,  because  the  purpose  of 
the  simulations  was  for  comparison  with  the  crushing  response  predicted  by  the  analytical  model  under 
specified  reflected  pressure  loading.  A  triangular  reflected  pressure  pulse  with  the  form  of  Eq.  (38a)  was 
applied  to  the  front  face  of  the  first  solid  element  in  the  computational  simulations.  As  shown  in  Table  1,  a 
nondimensional  peak  reflected  pressure  of  PR  =  10  was  used  in  both  the  “blast  pendulum”  and  “sandwich 
plate”  cases,  while  the  nondimensional  impulse  /R  was  selected  to  produce  comparable 
compaction  of  the  core  in  each  case,  and  a  significantly  larger  value  of  /R  was  used  in  the  “sandwich 
plate”  case. 

As  shown  in  Fig.  3,  the  computational  simulations  are  quite  sensitive  to  the  value  of  the  material  vis¬ 
cosity  coefficient  /(  specified  in  the  material  model.  For  the  “pendulum”  case  with  n  =  20  elements,  Fig.  3 
shows  that  the  computational  simulations  converge  very  closely  to  the  predictions  of  the  analytical  model 
when  /;  is  reduced  to  0.001  Pa  s,  while  larger  values  of  /i  result  in  reduced  core  deformation  and  increased 
values  of  the  nondimensional  back-face  acceleration  a2,  defined  in  Eq.  (31).  The  material  viscosity  coeffi¬ 
cient  /i  was  not  mentioned  by  Hanssen  et  al.  (2002),  but  that  study  did  find  that  more  than  400  elements 
were  required  to  achieve  adequate  convergence.  As  shown  in  Fig.  4,  for  a  given  value  of  /<  the  computa¬ 
tional  simulations  do  converge  toward  the  predictions  of  the  analytical  model  as  the  number  of  elements  is 
increased.  However,  with  /( =  0.05  Pa  s  (the  default  value  for  material  model  26),  the  computational  results 
have  not  yet  converged  even  with  n  =  150  elements,  which  is  consistent  with  the  observations  of  Hanssen 
et  al.  (2002).  By  reducing  p  to  0.001  Pa  s,  adequate  convergence  can  be  achieved  with  fewer  elements,  as 
shown  in  Fig.  3.  The  computational  results  in  Figs.  5  and  6  were  obtained  using  /(  =  0.001  Pas  and 
n=  150  elements,  which  was  found  to  be  more  than  sufficient  to  achieve  convergence  for  both  the  “blast 
pendulum”  and  the  “sandwich  plate”  cases  of  Table  1. 


Table  1 

Parameters  of  computational  simulations 


Case 

'7o 

h 

’ll 

Pa 

4 

Blast  pendulum 

0.0125 

0.0125 

0.975 

10 

0.015 

Sandwich  plate 

0.5 

0.25 

0.25 

10 

1 
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Fig.  3.  Influence  of  material  viscosity  on  convergence  of  computational  results  with  n  =  20  elements  (“blast  pendulum”  case),  (a) 
Nondimensional  core  thickness  and  (b)  nondimensional  back-face  acceleration.  Analytical  model,  o;  computational  results,  — . 


T  T 

Fig.  4.  Influence  of  number  of  elements  n  on  convergence  of  computational  results  with  /( =  0.05  Pa  s  (“blast  pendulum”  case),  (a) 
Nondimensional  core  thickness  and  (b)  nondimensional  back-face  acceleration.  Analytical  model,  o;  computational  results,  — . 


Fig.  5.  Time  histories  of  (a)  nondimensional  front-face  and  back-face  velocities  and  (b)  nondimensional  energy  quantities  for  “blast 
pendulum”  case.  Analytical  model,  o;  computational  results,  —  (n  =  150;  /i  =  0.001  Pa  s). 


In  Figs.  5  and  6,  computational  results  are  compared  with  predictions  of  the  analytical  model  for  the  “blast 
pendulum”  and  “sandwich  plate”  cases,  respectively,  and  in  both  cases,  good  agreement  is  observed  between 
the  analytical  and  computational  results.  Figs.  5a  and  6a  show  the  nondimensional  front-face  and  back-face 
velocities  defined  in  Eqs.  (28a, b),  and  Figs.  5b  and  6b  show  the  nondimensional  work  and  energy  quantities 
defined  in  Eqs.  (32a,b,c)  and  (36).  All  quantities  are  plotted  against  nondimensional  time  t,  defined  in  Eq.  (23), 


Fig.  6.  Time  histories  of  (a)  nondimensional  front-face  and  back-face  velocities  and  (b)  nondimensional  energy  quantities  for  “sandwich 
plate”  case.  Analytical  model,  o;  computational  results,  —  (n  =  150;  p  =  0.001  Pa  s). 


and  all  plots  terminate  at  the  instant  the  shock  front  is  arrested.  Beyond  this  instant,  the  entire  system  trans¬ 
lates  as  a  rigid  body.  At  the  instant  the  shock  wave  is  arrested,  the  velocities  of  the  front  and  back  faces  have 
converged,  so  that  tq  =  v2  =  1  for  both  the  “blast  pendulum”  and  the  “sandwich  plate”  cases,  and  the  non- 
dimensional  kinetic  energy  in  both  cases  has  plateaued  at  Ek  =  1.  While  the  final  value  of  the  nondimensional 
front-face  velocity  Si  is  the  same  for  both  the  “blast  pendulum”  and  the  “sandwich  plate”  cases,  comparison 
of  Figs.  5a  and  6a  shows  that  the  peak  value  of  v\  is  much  larger  in  the  “blast  pendulum”  case.  The  larger 
values  of  Vi  in  the  “blast  pendulum”  case  are  associated  with  much  larger  values  of  the  nondimensional  work 
done  by  the  pressure  pulse  Wp  than  in  the  “sandwich  plate”  case,  as  seen  by  comparison  of  Figs.  5b  and  6b. 
Figs.  5b  and  6b  also  show  that  the  nondimensional  energy  dissipation  Ed  is  significantly  larger  than  the  quasi¬ 
static  contribution  E{ j‘s  from  Eq.  (36),  demonstrating  that  the  rate-dependent  dissipation  associated  with  the 
shock  front  has  a  significant  influence  on  the  solution. 

4.  Initiation  and  arrest  of  densification 

In  order  to  initiate  densification  of  the  core,  the  peak  reflected  pressure  pK  must  be  sufficiently  high  that 
the  stress  behind  the  front  face  reaches  the  plateau  stress  erP.  The  maximum  pressure  than  can  be  sustained 
without  densification,  denoted  pp\  can  be  evaluated  by  considering  a  free-body  diagram  of  the  front  face 
under  the  applied  pressure  with  the  plateau  stress  erP  acting  at  the  interface  with  the  core.  If  densification 
has  not  initiated,  then  the  acceleration  of  the  front  face  must  equal  the  acceleration  of  the  center  of  mass, 
which  is  given  by  iiG  =  Pr'1  /m  according  to  Eq.  (1).  Newton’s  second  law  for  the  front  face  can  then  be  writ¬ 
ten  as  follows: 

pT  -  ffP  =  m^G  =  mipT/m)  (40) 

Eq.  (40)  can  be  rearranged  and  expressed  in  terms  of  nondimensional  variables,  yielding  the  following  expres¬ 
sion  for  the  nondimensional  peak  pressure  required  to  initiate  densification: 

=  (1  - 'll)'1  (41) 

If  Pr  >  Pjl".  then  densification  will  occur,  and  the  propagation  and  arrest  of  the  shock  front  can  be  evaluated 
by  solving  the  nondimensional  equation  of  motion  (25).  The  nondimensional  time  at  which  densification  is 
arrested,  denoted  rA,  can  be  evaluated  as  the  instant  that  the  nondimensional  shock-front  velocity  r  reaches 
zero.  The  value  of  the  nondimensional  coordinate  c  at  this  instant,  denoted  cA,  gives  the  mass  fraction  of 
uncrushed  core  material  remaining  at  arrest  of  densification. 

Figs.  7a  and  b  show  contours  of  qA  versus  PR  and  /R  for  the  “blast  pendulum”  and  “sandwich  panel”  cases, 
respectively  (see  Table  1),  with  a  triangular  pressure  pulse.  In  both  cases,  the  contours  for  £A  =  1  are  horizon¬ 
tal  lines  corresponding  to  Ppa  (41),  and  no  densification  occurs  for  PR  <  P'p".  The  contours  for  cA  =  0  in  Figs. 
7a  and  b  correspond  to  complete  densification  of  the  core.  For  values  of  PR  and  /R  beyond  these  contours,  the 
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Fig.  7.  Contours  of  remaining  mass  fraction  of  uncrushed 
nondimensional  reflected  pressure  PR  (triangular  pressure  pulse). 


nondimensional  impulse,  IR 


:  material  versus  nondimensional  reflected  impulse  IR  and 
"Blast  pendulum”  case  and  (b)  “sandwich  plate”  case. 


shock  front  is  not  arrested  before  reaching  the  back  face  (i.e.,  c  reaches  zero  before  u  reaches  zero),  and  the 
protective  function  of  the  cellular  core  is  lost,  as  discussed  previously.  The  minimum  nondimensional  impulse 
that  produces  complete  densification  of  the  core  is  denoted  7^nt,  and  the  £A  =  0  contours  in  Figs.  7a  and  b 
show  that  /))’"  increases  as  the  nondimensional  peak  pressure  PR  decreases. 

There  is  a  minimum  value  of  PK  for  which  complete  densification  can  be  achieved,  which  is  denoted  P'R\ 
When  PR  <  PR" ,  complete  densification  cannot  be  achieved,  even  if  the  pressure  is  maintained  indefinitely  (i.e., 
a  step  function,  for  which  IR  — >  oo).  Values  of  PRm  can  be  evaluated  numerically  by  selecting  a  very  large  value 
for  /R  (a  value  of  /R  =  106  was  found  to  be  sufficiently  large)  and  solving  Eq.  (25)  repeatedly  with  varying  PR 
to  find  the  minimum  value  of  PR  for  which  complete  densification  is  achieved.  In  the  “blast  pendulum”  case, 
P1™  =  1.019  while  P™*  =  1.013,  and  the  contours  for  cA  =  0  and  fA  =  I  in  Fig.  7a  are  almost  indistinguishable 
as  /R  becomes  large.  However,  in  the  “sandwich  plate”  case,  PRm  =  2.394  while  P™1  =  1.333,  and  a  clear  gap  is 
evident  between  the  contours  for  Ca  =  0  and  cA  =  I  as  /R  becomes  large. 

Plotted  in  grey  in  Figs.  7a  and  b  are  contours  corresponding  to  ta  =  rR,  where  rA  is  the  nondimensional 
time  at  which  the  shock  front  is  arrested  and  tr  is  the  nondimensional  duration  of  the  triangular  reflected  pres¬ 
sure  pulse.  For  values  of  PR  and  /R  below  the  rA  =  tr  contour,  the  shock  front  is  arrested  before  the  end  of  the 
pressure  pulse  (tA  <  tr),  and  the  response  is  more  sensitive  to  the  peak  pressure  than  to  the  total  impulse.  For 


Fraction  of  face-sheet  mass  in  front  face 
Hi  /  (h  +  %) 


Fraction  of  face-sheet  mass  in  front  face 
Hi  /(7i  +  Hi) 


Fig.  8.  Contours  with  varying  mass  distribution  of  (a)  nondimensional  reflected  pressure  required  to  initiate  densification  P™1  and  (b) 
minimum  nondimensional  reflected  pressure  for  which  complete  densification  of  the  core  can  be  achieved  P1™ . 
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values  of  FR  and  /R  above  the  ta  =  tr  contours,  the  shock  front  is  arrested  after  the  end  of  the  pressure  pulse 
(ta  >  tr),  and  the  response  is  more  sensitive  to  the  total  impulse  than  to  the  peak  pressure. 

The  dependence  of  F™‘  and  PR"  on  mass  distribution  is  shown  in  Figs.  8a  and  b,  respectively,  in  which  con¬ 
tours  of  F™'  and  Plj"'  are  plotted  against  the  fraction  of  mass  in  the  core  tj0,  and  against  the  fraction  of  remain¬ 
ing  mass  (i.e.,  face-sheet  mass)  in  the  front  face,  >]i/(>i\  +  i;2)-  The  remaining  figures  in  this  paper  use  the  same 
axes  as  Fig.  8  to  illustrate  the  influence  of  mass  distribution  on  various  quantities  of  interest.  To  aid  in  inter¬ 
preting  these  figures,  it  is  noted  that  the  “sandwich  plate”  case  corresponds  to  the  center  of  these  plots,  with 
half  of  the  mass  in  the  core,  i]0  =  0.5,  and  the  remaining  mass  distributed  equally  between  the  front  and  back 
faces,  >7i/(f/i  +  rj2)  =  0.5,  while  the  “blast  pendulum”  case  is  near  the  lower  left  corner  of  these  plots,  with 

=  0.0125  and  r\\l{r\\  + 1]2)  =  0.0127.  Fig.  8  shows  that  both  F™1  and  F1™  increase  with  increasing  rjJ 
O/i  +  '/eh  while  F™‘  decreases  with  increasing  )/0-  The  range  between  F™1  and  F^m,  over  which  an  infinite 
impulse  produces  only  partial  densification,  is  thus  maximized  in  the  upper  right-hand  corner  of  the  plots, 
which  corresponds  to  a  large  mass  fraction  in  the  core  with  most  of  the  remaining  mass  in  the  front  face. 

The  influence  of  mass  distribution  on  /£“,  the  nondimensional  reflected  impulse  required  for  complete  den¬ 
sification  of  the  core,  is  shown  in  Figs.  9a  and  b,  in  which  contours  of  /R“  for  FR  =  10  and  FR  =  100,  respec¬ 
tively,  are  plotted  using  the  same  axes  as  in  Fig.  8.  For  consistency  with  subsequent  calculations  that  include 
FSI  effects,  the  values  of  in  Fig.  9  were  computed  for  an  exponential  pressure  pulse  as  defined  in  Eq.  (39). 
It  was  previously  observed  in  Fig.  7  that  increases  with  decreasing  FR,  and  consequently,  larger  values  of 
7^nt  are  observed  for  FR  =  10  (Fig.  9a)  than  for  FR  =  100  (Fig.  9b).  While  the  contours  in  Figs.  9a  and  b  are 
fairly  close  over  much  of  the  plot,  significant  differences  are  observed  near  the  right-hand  edges  of  the  plots, 
which  correspond  to  a  large  fraction  of  the  face-sheet  mass  in  the  front  face.  In  Fig.  9a,  a  contour  is  plotted  for 
/))"  — >  oo,  which  is  the  same  as  the  contour  for  F^m  =  10  plotted  in  Fig.  8b.  Along  this  contour,  the  nondimen¬ 
sional  pressure  FR  =  10  must  be  sustained  indefinitely  in  order  to  achieve  complete  densification  of  the  core. 
The  shaded  region  beyond  this  contour  corresponds  to  F1™  >10,  and  in  this  region,  complete  densification  of 
the  core  cannot  be  achieved  with  FR  =  10,  even  for  an  infinite  impulse.  Fig.  9  shows  that  for  both  FR  =  10  and 
Fr  =  100,  7^nt  increases  with  increasing  +  /;2),  and  7^nt  generally  increases  with  increasing  rj 0,  except  in 
the  lower  right-hand  corner  of  Fig.  9a,  as  the  contour  for  /))’"  — >  oo  is  approached. 

Figs.  10a  and  b  show  the  influence  of  mass  distribution  on  tA,  the  nondimensional  time  at  which  the  shock 
front  is  arrested,  for  an  exponential  reflected  pressure  pulse  with  FR  =  10  and  FR  =  100,  respectively.  Con¬ 
tours  of  ta/tr  are  plotted,  where  tr  =  1/FR  is  the  nondimensional  decay  period  of  the  exponential  pulse. 
To  enable  a  consistent  comparison,  since  ta  depends  on  /R,  the  values  in  Figs.  10a  and  b  were  evaluated  using 
the  minimum  value  of  /R  for  which  complete  densification  is  achieved  in  each  case  (i.e.,  the  values  of  7™* 
shown  in  Figs.  9a  and  b,  respectively).  In  the  shaded  region  of  Fig.  9b,  for  which  complete  densification  cannot 
be  achieved,  a  value  of  7R  =  100  was  used.  Fig.  10  shows  that  the  values  of  ta/tr  are  about  10  times  larger  for 


Fraction  of  face-sheet  mass  in  front  face  Fraction  of  face-sheet  mass  in  front  face 

7 1  I  (?7i  +  Vi)  i]l  /  (?/,  +  r/2) 

Fig.  9.  Contours  with  varying  mass  distribution  of  nondimensional  reflected  impulse  /)(“  required  for  complete  densification  of  core 
(exponential  pressure  pulse),  (a)  PR  =  10  and  (b)  PR  =  100. 


J.A.  Main,  G.A.  Gazonas  /  International  Journal  of  Solids  and  Structures  45  (2008)  2297—2321 


2311 


Fraction  of  face-sheet  mass  in  front  face  Fraction  of  face-sheet  mass  in  front  face 

%/(.%  +  %)  Hi  l(Th  +  %) 

Fig.  10.  Contours  with  varying  mass  distribution  of  the  ratio  ta/tr  between  the  time  of  arrest  of  densification  and  the  decay  period  of  the 
reflected  pressure  pulse  (values  correspond  to  complete  densification  of  the  core,  with  reflected  impulses  shown  in  Fig.  9).  (a)  PR  =  10  and 
(b)  PR  =  100. 

PR  =  100  than  for  PR  =  10,  as  a  consequence  of  the  smaller  value  of  tr  in  the  former  case.  For  both  PR  =  10 
and  PR  =  100,  tA/rR  decreases  with  increasing  t;0  and  with  increasing  r\J(r]\  +  jj2).  Fig.  10a  shows  that  in  some 
cases  ta  <  tr.  It  was  observed  previously  in  Fig.  7  that  when  tA  <  tr,  the  response  becomes  more  sensitive  to 
the  peak  pressure  than  to  the  total  impulse.  From  Fig.  10a  it  can  then  be  observed  that  when  t j0  and/or  t; j/ 
('h  +  Ui)  become  large  (i.e.,  beyond  the  contour  for  tA  <  tR  =  1),  the  system  becomes  more  sensitive  to  the 
peak  pressure  and  less  sensitive  to  the  total  impulse.  In  Fig.  10b  a  contour  for  ta/tr  =  0  is  shown,  which  is 
the  same  as  the  contour  for  75™1  =  10  shown  previously  in  Fig.  8a.  In  the  shaded  area  below  this  contour, 
the  nondimensional  peak  pressure  of  PR  =  10  is  insufficient  to  initiate  densification,  regardless  of  the  applied 
impulse. 

Fig.  1 1  shows  contours  with  varying  mass  distribution  of  the  ratio  EA/EAS  of  the  total  energy  dissipated  in 
the  core  to  the  quasi-static  energy  dissipation  given  by  Eq.  (36).  As  in  Fig.  10,  the  contours  in  Fig.  11  corre¬ 
spond  to  complete  densification  of  the  core  and  were  computed  using  the  values  of  7^nt  shown  in  Fig.  9,  with 
7R  =  100  for  cases  in  which  complete  densification  cannot  be  achieved.  Fig.  11  shows  that  Ea/E^s  increases 
with  increasing  tj0  for  both  PR  =  10  and  PR  =  100,  while  significantly  larger  values  of  Ea/EY  are  observed 


7i  /  (7i  +  Ih)  m  /  (7i  +  'll) 


Fig.  11.  Contours  with  varying  mass  distribution  of  the  ratio  Ea/EY  between  the  total  energy  dissipated  in  the  core  and  the  quasi-static 
energy  dissipation  (values  correspond  to  complete  densification  of  the  core,  with  reflected  impulses  shown  in  Fig.  9).  (a)  PR  =  10  and  (b) 
PR=  100. 
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Fig.  12.  Influence  of  mass  distribution  for  a  Dirac  delta  pulse  (PR  — »  oo):  contours  of  (a)  nondimensional  reflected  impulse  /jj[“  required 
for  complete  densification  of  core  and  (b)  ratio  EA/E^  between  the  total  energy  dissipated  in  the  core  and  the  quasi-static  energy 
dissipation  (values  correspond  to  complete  densification  of  the  core,  with  reflected  impulses  shown  in  (a)). 


for  PR  =  100.  For  large  ;/0,  the  rate-dependent  contribution  to  the  energy  dissipation  can  be  many  times  larger 
than  the  quasi-static  contribution.  The  shaded  region  in  the  lower  right-hand  corner  of  Fig.  11a  corresponds 
to  the  shaded  region  in  Fig.  10a,  in  which  the  nondimensional  peak  pressure  of  Fr  =  10  is  insufficient  to  ini¬ 
tiate  densification  of  the  core,  regardless  of  the  applied  impulse.  No  energy  is  dissipated  in  this  region,  so  the 
ratio  Ea/EY  is  undefined. 

While  Figs.  9-11  compared  results  for  two  different  levels  of  peak  pressure  (FR  =  10  and  PR  =  100), 
Fig.  12  presents  results  for  the  limiting  case  of  a  Dirac  delta  pulse,  for  which  PR  — >  oo  and  iR  — >  0.  In  this 
limiting  case,  the  values  of  nondimensional  peak  pressure  required  to  initiate  and  complete  densification 
become  irrelevant  [P'Rl  and  shown  in  Fig.  8  are  always  exceeded),  and  complete  densification  can  always 
be  achieved  for  a  finite  value  of  the  nondimensional  reflected  impulse  /R.  Fig.  12a  presents  contours  with 
varying  mass  distribution  of  the  critical  nondimensional  impulse  /Rm  required  for  complete  densification. 
Note  that  the  contours  of  /Rr,t  for  the  Dirac  delta  pulse  in  Fig.  12a  are  virtually  indistinguishable  from  those 
for  PR  =  100  plotted  previously  in  Fig.  9b.  Although  not  plotted  herein,  it  is  noted  that  values  of  the  non- 
dimensional  time  ta  required  to  arrest  densification  for  the  Dirac  delta  pulse  are  also  nearly  equivalent  to 
those  for  PR  =  100.  (The  ratio  iA/tR,  plotted  for  /JR  =  10  and  PR  =  100  in  Fig.  10,  is  undefined  for  the  Dirac 
delta  pulse  because  tr  — >  0.)  The  energy  dissipation  in  the  core,  however,  can  be  significantly  larger  for  the 
Dirac  delta  pulse  than  for  PR  =  100,  as  can  be  observed  by  comparing  the  contours  of  Ea/EY  for  the  Dirac 
delta  pulse  in  Fig.  12b  with  those  shown  previously  for  PR  =  100  in  Fig.  lib.  When  the  mass  fraction  in  the 
core  1(0  approaches  unity,  the  rate-dependent  contribution  to  the  energy  dissipation  can  be  significantly  lar¬ 
ger  for  the  Dirac  delta  pulse  than  for  PR  =  100.  These  observations  are  consistent  with  those  of  Vaziri  and 
Hutchinson  (2006)  in  computational  simulations  of  blast  loading  on  metal  sandwich  plates,  who  noted  sig¬ 
nificantly  larger  energy  dissipation  when  the  blast  impulse  was  imparted  as  an  initial  velocity  to  the  front 
face  than  when  it  was  applied  as  an  exponential  pressure  pulse,  even  though  the  back-face  deflections  in 
the  two  cases  were  fairly  close. 


5.  Fluid-structure  interaction 

The  analytical  model  in  Section  2  was  formulated  in  terms  of  the  reflected  pressure  pT{t),  which  is  influenced 
by  FSI  effects.  In  this  section,  the  approximation  developed  by  Kambouchev  et  al.  (2006)  is  used  to  define  the 
reflected  pressure  pft)  corresponding  to  a  specified  incident  pressure  pulse  pft),  accounting  for  nonlinear  com¬ 
pressibility  and  FSI  effects.  For  consistency  with  Kambouchev  et  al.  (2006),  an  exponential  incident  pressure 
pulse  is  considered,  of  the  form  p{(t)  =  p0e~(/'°,  where  t0  is  the  decay  period.  Kambouchev  et  al.  (2006) 
observed  that  the  reflected  pressure  pft)  can  also  be  captured  quite  well  by  an  exponential  pulse,  which  is  con- 
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veniently  expressed  as  pr{t)  =  /?Re  '/<r,  where  ?R  is  the  decay  period  of  the  reflected  pulse.  The  peak  reflected 
pressure  pR  is  related  to  the  peak  incident  pressure  through  the  pressure  reflection  coefficient  CR: 


_Pr  =  27  +  4(Po/Pa) 
-  Po  1  +  Oo/Ra) 


(42) 


in  which  pA  is  the  ambient  air  pressure. 

The  key  result  of  Kambouchev  et  al.  (2006)  is  an  approximation  relating  the  reflected  impulse  zR  to  the  inci¬ 
dent  impulse  i0,  which  can  be  expressed  as  an  impulse  reflection  coefficient  aR: 


*r  =  —  ^  7r 
io 


&/0+&) 


(43) 


where  yR  and  /R  are  defined  as  follows: 


yR  =  8-42^1n(l+^)  (44) 

f  =  (aEi+i\  I  (6  +  Cr)(p0/^a)  +  7  (  , 

■  V  Pa  Jy(Po/PA  +  mi  +  6C*)(p0/pA)  +  T][C*(p0/pA)  +  T]  [  1 


The  nondimensional  FS1  parameter  /?s  in  Eq.  (43)  is  conveniently  expressed  as 


Ps 


Pa 


n+6(po/pA)\ 
V  1+Po/Pa  ) 


(46) 


where  /1a  is  an  alternative  nondimensional  parameter  introduced  by  Vaziri  and  Hutchinson  (2006),  expressed 
in  terms  of  ambient  air  properties: 

Pa  =  (47) 

wfsi 

in  which  pA  is  the  ambient  air  density,  cA  =  a/1.4 pA/pA  is  the  ambient  speed  of  sound  in  air,  and  /nFSi  is  the 
areal  density  of  the  plate  considered  in  the  FSI  analysis. 

Kambouchev  et  al.  (2006)  considered  a  freestanding  solid  plate  in  developing  Eq.  (43),  and  further  approx¬ 
imations  must  be  introduced  in  selecting  an  appropriate  value  of  z;7FSI  to  extend  this  result  to  a  sandwich  plate 
with  a  crushable  core.  Vaziri  and  Hutchinson  (2006)  considered  only  the  mass  of  the  front  face  in  their  analysis 
of  FSI  for  sandwich  plates,  i.e.,  setting  mFSI  =  m\  in  Eq.  (47).  However,  they  noted  that  this  neglects  the  resis¬ 
tance  due  to  crushing  of  the  core  and  thus  overestimates  the  beneficial  effects  of  FSI.  In  this  paper,  an  alter¬ 
native  approximation  is  proposed  that  incorporates  the  resistance  of  the  core  by  defining  ;«FSi  as  the  areal 
density  of  a  solid  plate  with  the  same  initial  acceleration  of  the  front  face  of  the  sandwich  plate. 

Provided  that  the  peak  reflected  pressure  pR  is  sufficiently  high  to  initiate  densification,  the  initial  acceler¬ 
ation  of  the  front  face  of  the  sandwich  plate  is  given  by  m™1  =  {pR  -  erP)/w|.  Requiring  a  corresponding  solid 
plate  to  have  the  same  initial  acceleration  under  the  applied  pressure  pR,  mFSi  can  then  be  defined  as 
'«fsi  =  Pr/'^T1  =  "h/(l  ~  o'p/z’r).  ^  Pr  is  insufficient  to  initiate  densification,  then  the  sandwich  plate 
responds  as  a  monolithic  solid  plate  and  mFsi  =  m.  These  results  can  be  expressed  in  terms  of  nondimensional 
variables  as 


_  mFSl 

Z/FSI  =  -  — 

m 


>7i(l  ~PrT\  Pr>PT 
1,  Pr<PT 


(48) 


where  )fFS i  represents  the  effective  mass  fraction  of  the  sandwich  panel  that  participates  in  FSI,  and  PRli  is  de¬ 
fined  in  Eq.  (41).  Eq.  (48)  shows  that  z/FSi  — >  ip  as  the  nondimensional  peak  pressure  Pr  becomes  large,  which 
corresponds  to  the  approximation  used  by  Vaziri  and  Hutchinson  (2006). 

To  check  the  appropriateness  of  using  Eq.  (48)  to  extend  FSI  results  for  solid  plates  to  sandwich  plates,  the 
acceleration  of  the  solid  plate  considered  in  the  FSI  analysis,  denoted  i/FSI  =  pr(t) /wFSI,  can  be  compared  with 
the  acceleration  of  the  front  face  of  the  sandwich  panel,  given  by  ii\  =  \pr(t)  —  ffp]/[wi  +  p0(4  —  *)].  With  the 
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exponential  reflected  pressure  pulse p,.(t)  =  pRe  and  with  mFSI  defined  by  Eq.  (48),  the  ratio  mfsi/mi  can  be 
expressed  in  terms  of  nondimensional  variables  as  a  product  of  two  factors: 

^FSI  _  >/l  +  >7p(  1  ~  0  _  1  ~  ^R1 

hi  m  1  -  Pr  eT/TR  ^ 

By  definition,  hFS i/u\  =  1  initially,  but  both  factors  in  Eq.  (49)  increase  as  time  progresses  and  c  decreases 
from  unity.  The  first  factor  in  Eq.  (49)  is  associated  with  the  accretion  of  crushed  core  material  as  the  densi- 
fication  front  propagates.  The  resulting  increase  in  inertia  may  significantly  reduce  the  front-face  accelerations 
if  the  mass  of  the  front  face  is  small  relative  to  the  mass  of  the  core.  The  second  factor  in  Eq.  (49)  is  associated 
with  the  resistance  provided  by  the  plateau  stress  of  the  core  material,  which  remains  constant  as  the  reflected 
pressure  decays.  The  ratio  hFSI/«i  is  undefined  for  t  /tr  =  In  /JR,  because  the  plateau  stress  equals  the  reflected 
pressure  at  this  instant,  and  therefore  hi  =  0.  For  t  /tr  >  In  PR,  the  plateau  stress  exceeds  the  reflected  pres¬ 
sure,  and  the  front  face  decelerates  (hi  <  0),  provided  that  the  densification  front  has  not  already  been  ar¬ 
rested.  Both  of  these  effects  (the  accretion  of  crushed  core  material  and  the  resistance  provided  by  the 
plateau  stress)  lead  to  reduced  accelerations  of  the  front  face  relative  to  the  accelerations  of  the  solid  plate 
considered  in  the  FSI  analysis.  If  such  reductions  are  significant  during  the  period  of  application  of  the  pres¬ 
sure  pulse,  then  Eq.  (48)  will  overestimate  the  motion  of  the  front  face  and  thus  overestimate  the  beneficial 
effects  of  FSI.  Conversely,  if  the  ratio  hFSi/hi  remains  fairly  close  to  unity  during  the  application  of  the  pres¬ 
sure  pulse,  then  Eq.  (48)  should  be  reasonably  accurate. 

Figs.  13a  and  b  show  contours  with  varying  mass  distribution  of  the  ratio  hFSI(f*)/iii(c)for  PR  =  10  and 
PR  =  100,  respectively.  The  time  f*  is  the  lesser  of  the  decay  period  fR  and  the  arrest  time  tA.  so  that  the  ratio 
hFS i/iii  is  evaluated  at  t  =  tR  unless  the  densification  front  is  arrested  first  (see  the  contours  of  tA/tR  in  Fig.  10). 
Note  that  at  t  =  tR,  the  reflected  pressure  has  decayed  to  36.8%  of  its  initial  value,  and  63.2  %  of  the  total 
impulse  has  been  imparted.  As  in  Figs.  10  and  11,  the  contours  in  Fig.  13  correspond  to  complete  densification 
of  the  core  and  were  computed  using  the  values  of  IRl  shown  in  Fig.  9,  with  /R  =  100  for  cases  in  which  com¬ 
plete  densification  cannot  be  achieved.  When  the  mass  fraction  in  the  core  rj0  is  small,  the  ratio  «FSi/wiis  fairly 
close  to  unity  for  both  PR  =  10  and  PR  =  100.  However,  large  values  of  iiFSi/iii  are  observed  as  t]0  approaches 
unity,  particularly  when  the  mass  fraction  in  the  front  face  is  small.  Much  larger  values  of  iiFS i/u\  are  observed 
for  PR  =  10  than  for  PR  =  100,  indicating  that  it  is  less  appropriate  in  the  former  case  to  approximate  FSI 
effects  using  Eq.  (48).  This  is  due  largely  to  the  longer  decay  period  associated  with  the  lower  level  of  peak 
pressure,  which  allows  more  time  for  the  accretion  of  crushed  core  material  during  the  application  of  the  pres¬ 
sure  pulse.  For  this  reason,  in  the  subsequent  analysis  of  FSI  effects,  numerical  results  are  presented  only  for 
the  higher  level  of  peak  pressure,  PR  =  100.  In  addition,  the  contour  for  iiFsi(f*)/“i(**)  =  1.5  in  Fig.  13b  is 


%  I  Vh  +  >h)  Vi  /  ill  +  %) 


Fig.  13.  Contours  with  varying  mass  distribution  of  the  ratio  i/Fsi(t»)/ui(n)  between  the  acceleration  of  the  “equivalent”  solid  plate  used 
in  treatment  of  FSI  and  the  actual  front-face  acceleration  of  the  sandwich  plate.  The  time  t,  is  the  lesser  of  the  decay  period  tR  and  the 
arrest  time  tA  (see  Fig.  10).  (a)  PR  =  10  and  (b)  PR  =  100. 
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tentatively  identified  as  an  approximate  limit  of  applicability  of  Eq.  (48),  and  in  subsequent  plots  in  which  FSI 
effects  are  considered,  the  region  above  this  contour  is  hatched  to  indicate  that  the  approximation  in  Eq.  (48)  is 
not  appropriate.  Note  that  this  region  corresponds  to  small  values  of  rj  i  and/or  large  values  of  tj0,  so  that  the 
mass  of  the  core  is  large  relative  to  the  mass  of  the  front  face. 

Implementing  the  FSI  approximations  in  this  section  within  the  analytical  model  of  Section  2  first  requires 
specification  of  a  nondimensional  incident  pressure  P0  and  a  nondimensional  incident  impulse  70,  defined  con¬ 
sistently  with  the  corresponding  reflected  quantities  in  Eqs.  (37)  and  (26): 


Po  = 


Po_. 

1 

<Tp 


*0  /  Po 

m  Y  ffp£D 


(50a,  b) 


The  FSI  approximations  depend  on  po/pA,  which  can  be  expressed  in  terms  of  Po  as  p<Jp,\  =  {nflp.flPo-  In  the 
numerical  results  presented  subsequently,  a  value  of  r,flpA  =  10  is  considered,  which  corresponds  to  a  plateau 
stress  of  about  erP  =  1  MPa,  as  considered  previously  in  the  computational  simulations.  With  po/p,\  thus  spec¬ 
ified,  the  pressure  reflection  coefficient  CR  (42)  can  be  evaluated,  and  the  nondimensional  peak  reflected  pres¬ 
sure  can  be  computed  as  PR  —  CRP0. 

It  follows  from  Eqs.  (50a, b)  that/o/Po  =  (to/»?)\/po<Tp/eD,  whereby  /?A  (47)  can  be  expressed  in  terms  of  P0 
and  70  as 

_  PaCa  1  7q  (-51) 

JpJrpfz d  *1 fsi  p o 

where  ^FSi  can  be  evaluated  from  Eq.  (48),  and  a  value  of  the  nondimensional  constant  pAcA/  Jpo^pJ^d  must 
be  specified  from  the  relevant  physical  properties.  In  the  numerical  results  presented  subsequently,  a  value  of 
Paca/v/Po<tp/£d  =  0.02206  is  considered,  which  corresponds  to  setting  pA  =  1.225  kg/nr3  and  cA  =  340.3  m/s 
for  air  and  using  the  values  of  p0  =  250  kg/m3,  erP  =  1  MPa,  and  e0  =  0.7  considered  previously  in  the  compu¬ 
tational  simulations.  The  nondimensional  FSI  parameter  /is  (46)  and  the  impulse  reflection  coefficient  aR  (43) 
can  then  be  evaluated  in  terms  of  p0/pA  and  [SA,  whereby  the  nondimensional  reflected  impulse  can  be  obtained 
as  IR  =  aR70.  With  PR  and  IR  thus  determined,  the  response  of  the  sandwich  plate  can  be  evaluated  by  solving 
the  nondimensional  equation  of  motion  (25)  with  the  exponential  nondimensional  reflected  pressure  pulse  in 
Eq.  (39). 

Fig.  14a  shows  contours  with  varying  mass  distribution  of  7jjnt,  the  nondimensional  incident  impulse 
required  to  produce  complete  compaction  of  the  core  for  PR  =  100.  The  values  of  in  Fig.  14a  are  related 
to  the  values  of  in  Fig.  9b  through  the  impulse  reflection  coefficient  aR  (43):  7™*  =  ocrZ™1.  Fig.  14b  shows 


Fraction  of  face-sheet  mass  in  front  face 
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Fig.  14.  Contours  with  varying  mass  distribution  of  (a)  the  nondimensional  incident  impulse  /™‘  required  for  complete  densification  of  the 
core  and  (b)  the  corresponding  ratio  SrA/r  between  the  impulse  reflected  from  the  sandwich  plate  and  the  impulse  reflected  from  a  fixed, 
rigid  wall.  The  approximate  treatment  of  FSI  is  not  applicable  in  the  hatched  region,  in  which  iiFSi(t»)/ifi(t.)  >1.5  (_PR  =  100,  P0  =  13 ,p0/ 
pA  =  130,  pAcA/y/p0aP/ED  =  0.02206). 
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corresponding  contours  of  the  ratio  aR/yR  between  the  impulse  reflected  from  the  sandwich  plate  and  the 
impulse  reflected  from  a  fixed,  rigid  wall  for  PK  =  100.  The  beneficial  effects  of  FSI  are  evident  in  Fig.  14b, 
with  values  of  aR/yR  less  than  unity  indicating  impulse  reductions  due  to  FSI.  Fig.  14b  shows  a  trend  of 
decreasing  aR/yR  with  increasing  rj0,  with  impulse  reductions  of  more  than  30  %  observed  within  the  domain 
of  applicability  of  the  FSI  approximation  (48).  It  is  noted  that  these  decreases  in  aR/yR  are  associated  with 
underlying  increases  in  the  nondimensional  FSI  parameter  />A  (43).  Reductions  of  aR/yR  with  increasing  /iA 
follow  from  Eq.  (43)  and  are  illustrated  in  Fig.  1  of  Vaziri  and  Hutchinson  (2006).  (The  ratio  aR/yR  is  equiv¬ 
alent  to  I/Im p-,,,0  in  their  notation.) 

The  increases  in  /1A  that  underlie  Fig.  14b  result  from  two  effects.  Firstly,  the  contours  in  Fig.  14b  corre¬ 
spond  to  complete  densification  of  the  core,  with  the  critical  values  of  incident  impulse  7™'  shown  in  Fig.  14a. 
Increases  in/™'  then  produce  increases  in  f>A.  because  /iA  is  proportional  to  I0  as  shown  by  Eq.  (51).  Secondly, 
decreases  in  produce  increases  in  /?A,  because  /?A  is  inversely  proportional  to  i]u  according  to  Eqs.  (51)  and 
(48).  This  latter  effect  leads  to  reductions  in  aR/yR  as  the  fraction  of  face-sheet  mass  in  the  front  face  r\\! 
Oh  +  >12)  decreases,  due  to  the  relatively  smaller  mass  participating  in  FSI.  It  is  a  consequence  of  this  latter 
effect  that  the  contours  of  7™'  in  Fig.  14a  curve  downwards  as  they  approach  the  hatched  region  along  the 
left-hand  edge  of  the  plot,  in  contrast  with  the  corresponding  contours  of  7^"  in  Fig.  9b.  Further  increases 
in  7™'  may  occur  in  hatched  region  where  iq /( )/ 1  +  i/2)  becomes  small.  However,  in  the  domain  of  applicability 
of  the  FSI  approximation  (48),  7™'  increases  with  increasing  ;/i/()/]  +  >72),  and  the  contours  for  7™'  follow  the 
same  general  trend  as  those  for  7^“  in  Fig.  9b. 

It  is  interesting  to  compare  the  work  done  by  the  pressure  pulse  with  the  work  done  on  a  corresponding 
solid  plate  with  the  same  total  mass.  For  a  prescribed  reflected  pressure  pulse  pr(  t),  Eq.  (12)  shows  that  crush¬ 
ing  of  the  cellular  core  enhances  the  work  done  relative  to  that  done  on  a  solid  plate  due  to  the  increased  veloc¬ 
ity  of  the  front  face.  Due  to  the  beneficial  effects  of  FSI,  however,  pT(t)  is  not  simply  prescribed,  but  the 
reflected  impulse  for  a  sandwich  plate  can  be  reduced  by  the  increase  motion  of  the  front  face.  Combining 
Eqs.  (32a)  and  (43),  the  total  work  done  by  the  pressure  pulse  on  the  sandwich  plate  as  a  consequence  of  these 
competing  effects  can  be  expressed  as  fFp0t  =  !Fp0t^  (aRi0)2/m\,  where  lT‘pot  is  evaluated  by  integrating  Eq.  (34), 
and  the  superscript  “tot”  denotes  limiting  values  as  t  — >  00.  The  work  done  on  a  solid  plate  with  the  same  total 
mass  can  be  expressed  as  IFp0'olid  =  { (aRlKii0)2/m,  where  a^lld  is  evaluated  from  Eq.  (43)  with  mFSI  =  m  (or 
f/Fsr  =1).  The  ratio  of  the  work  done  on  the  sandwich  plate  to  the  work  done  on  the  corresponding  solid  plate 
can  then  be  expressed  as 


W 

_ p _ =  If tot 

M/ t°t  P 

vy  p, solid 


(52) 


Fraction  of  face-sheet  mass  in  front  face 
Hi  /  (  m  +  >k) 


Fig.  15.  Contours  with  varying  mass  distribution  of  the  ratio  ^p'/^pLud  between  the  work  done  by  the  pressure  pulse  on  the  sandwich 
plate  and  the  work  done  on  a  solid  plate  with  the  same  total  mass.  Values  correspond  to  complete  densification  of  the  core,  with  incident 
impulses  shown  in  Fig.  14a.  The  approximate  treatment  of  FSI  is  not  applicable  in  the  hatched  region,  in  which  «Fsi(t*)/“i(f.)  >1.5 


(PR  =  100,  P0  =  13,  Po/pa  =  130,  Pa'WvVo'U'Ad  =  0.02206). 
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Fig.  15  shows  contours  with  varying  mass  distribution  of  the  ratio  Wl°l  /  Wl°\oXiA  f°r  =  100.  As  in  Fig.  14b, 
the  contours  in  Fig.  15  correspond  to  complete  densification  of  the  core,  with  the  critical  values  of  incident 
impulse  7™*  shown  in  Figs.  14a.  All  of  the  values  in  Fig.  15  are  greater  than  unity,  indicating  that  impulse 
reductions  due  to  the  beneficial  effects  of  FSI  are  insufficient  to  offset  the  enhanced  work  that  results  from 
increased  velocities  of  the  front  face.  Fig.  15  shows  that  /  Wl°\oXiA  is  largest  in  the  lower  left  corner  of 
the  plot,  corresponding  to  a  small  mass  fraction  in  the  core  with  most  of  the  mass  in  the  back  face.  Interest¬ 
ingly,  this  mass  distribution  corresponds  to  the  blast  pendulum  experiments  discussed  previously,  for  which 
simulation  results  were  shown  in  Fig.  5.  Fig.  15  shows  that  for  such  mass  distributions,  the  work  done  by 
the  pressure  pulse  can  be  more  than  ten  times  larger  than  the  work  done  on  a  corresponding  solid  plate. 
Increasing  the  fraction  of  face-sheet  mass  in  the  front  face  171/(1/1  +  i/2)  leads  to  decreases  in  Wl°l /  Wl°loXiA. 

6.  Mitigation  of  back-face  accelerations 

The  nondimensional  back-face  acceleration  a2,  defined  in  Eq.  (31),  is  an  indication  of  the  degree  of  protec¬ 
tion  provided  by  the  limiting  plateau  stress  uP  of  the  cellular  core.  According  to  Eq.  (31),  the  peak  back-face 
acceleration  occurs  when  <f  =  0  (i.e.,  at  the  instant  the  shock  front  reaches  the  back  face)  and  is  given  by 
a™*  =  fjf '  ’  while  the  minimum  back-face  acceleration  occurs  when  c  =  1  (i.e.,  at  the  instant  that  densification 
initiates)  and  is  given  by  a™1  =  (tj0  +  rj2)~l.  Figs.  16a  and  b  show  contours  of  a'2'1  and  S™1,  respectively,  with 
varying  mass  distribution.  The  two  quantities  are  equivalent  in  the  limit  as  i/0  — >  0,  but  as  increases, 
decreases  while  a™1  increases.  For  a  given  mass  fraction  in  the  core  i/0,  increasing  the  mass  fraction  in  the  front 
face  leads  to  increases  in  both  ei'”"  and  a2'\  thus  producing  larger  back-face  accelerations. 

While  the  nondimensional  back-face  acceleration  a2  (31)  is  defined  in  terms  of  the  plateau  stress  <xP,  an 
alternative  nondimensional  back-face  acceleration  can  be  defined  as  ii2/ u oak,  where  =  pR/m  is  the  peak 

acceleration  of  the  center  of  mass.  (This  alternative  nondimensionalization  cannot  be  used  for  the  Dirac  delta 
pulse,  for  which  pR  — >  oo  and  thus  is  undefined.)  It  follows  from  Eq.  (31)  that  u2/u^:ik  =  a2/PR,  so  that 
values  of  iii/u^A  can  be  obtained  from  corresponding  values  of  a2  simply  through  division  by  the  nondimen¬ 
sional  peak  reflected  pressure  PK.  Using  this  nondimensionalization,  the  peak  back-face  acceleration  (at  com¬ 
plete  densification)  can  be  expressed  as  u2'l/iu^ik  =  a^/P^  =  (r]2P^)  \  Interestingly,  when  rj2  <  UPr  (or 
equivalently,  when  a™4  >  Pr),  the  peak  back-face  acceleration  is  actually  larger  than  the  corresponding  peak 
rigid-body  acceleration,  thus  sacrificing  the  protective  role  of  the  cellular  core  in  limiting  the  back-face  face 
accelerations.  With  7^=10,  for  example,  the  contour  for  a2A  =  50  in  Fig.  12b  corresponds  to 
u2"-/u^A  =  5,  meaning  that  the  peak  back-face  acceleration  is  five  times  larger  than  the  corresponding 
rigid-body  acceleration. 


mKm+TJi)  ihKrh  +  Th) 

Fig.  16.  Contours  with  varying  mass  distribution  of  nondimensional  back-face  acceleration,  a2.  (a)  At  initiation  of  densification,  and 
(b)  at  complete  densification  of  core,  o™‘. 
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The  plots  shown  previously  in  Figs.  8b,  9,  14a,  and  15  indicate  that  for  a  given  mass  fraction  in  the  core  r\Q, 
increasing  the  mass  fraction  in  the  front  face  increases  the  resilience  of  the  system  by  increasing  the  nondimen- 
sional  peak  pressure  P Rm  and  the  nondimensional  reflected  and  incident  impulses,  /Rm  and  7™',  required  for 
complete  densification  of  the  core,  and  also  by  reducing  the  nondimensional  work  /W^soM  done  by  the 
pressure  pulse  relative  to  that  done  on  a  corresponding  solid  plate  with  the  same  total  mass.  However, 
Fig.  16  indicates  that  increasing  the  mass  fraction  in  the  front  face  also  leads  to  increased  back-face  acceler¬ 
ations,  thus  sacrificing  a  protective  function  of  the  cellular  core.  A  design  optimization  problem  can  then  be 
posed  by  seeking  to  determine  the  distribution  of  mass  that  maximizes  either  the  nondimensional  reflected 
impulse  7R  or  the  nondimensional  incident  impulse  70  that  can  be  sustained,  while  limiting  the  back-face  accel¬ 
erations  to  a  specified  maximum  value,  denoted  5“ax. 

Figs.  17a  and  b  show  contour  plots  with  varying  mass  distribution  of,  respectively,  the  maximum  nondi¬ 
mensional  reflected  impulse  7R  and  the  maximum  nondimensional  incident  impulse  70  that  can  be  sustained 
with  a  maximum  allowable  back-face  acceleration  of  a”ax  =  5.  The  values  of  70  in  Fig.  17b  are  related  to 
the  values  of  /R  in  Fig.  17a  through  the  impulse  reflection  coefficient  aR  (43),  while  values  of  70  are  not  pre¬ 
sented  in  the  hatched  region  of  Fig.  17b  because  the  FS1  approximation  (48)  is  not  appropriate.  The  plots  in 
Figs.  17a  and  b  can  each  be  divided  into  three  regions  in  which  different  conditions  prevail.  The  grey  curve  in 
each  figure  corresponds  to  the  contour  for  ac2nt  =  5,  plotted  previously  in  Fig.  16b.  In  the  region  below  the  grey 
curve  ii2"  <  5,  so  complete  densification  of  the  core  can  be  permitted  without  exceeding  5™ax.  Therefore,  in  this 
region  the  values  of  maximum  permissible  impulse  correspond  to  complete  densification  of  the  core,  and  the 
contours  in  Fig.  17a  are  the  same  as  those  for  7 in  Fig.  9b,  while  the  contours  in  Fig.  17b  are  the  same  as 
those  for  7™‘  in  Fig.  14a.  In  the  region  above  the  grey  curve  a™1  >  5,  so  a™ax  will  be  exceeded  if  complete  den¬ 
sification  is  allowed.  Therefore,  the  maximum  permissible  impulses  in  this  region  correspond  to  partial  den¬ 
sification  of  the  core  and  the  values  in  Figs.  17a  and  b  are  less  than  7™'  and  7™‘,  respectively.  The  shaded 
region  in  the  lower  right-hand  corner  of  Figs.  17a  and  b  is  bounded  above  by  the  contour  for  =  5,  plotted 
previously  in  Fig.  16a.  In  this  region  a™1  >  5,  so  a”ax  will  be  exceeded  if  densification  is  even  allowed  to  ini¬ 
tiate,  and  thus  a  nonzero  impulse  cannot  be  permitted. 

Figs.  17a  and  b  show  that  the  allowable  reflected  impulse  7R  and  the  allowable  incident  impulse  70  both 
increase  monotonically  with  increasing  mass  fraction  in  the  core  i/0.  It  can  also  be  observed  in  Figs.  17a 
and  b  that  for  a  given  )/0,  the  allowable  values  of  7R  and  70  are  maximized  along  the  contour  for  a™1  =  5, 
shown  in  grey.  This  is  evident  from  the  fact  that  the  contours  of  7R  and  70  attain  their  minimum  values  of 
/; o  along  the  grey  curve,  with  attendant  changes  in  sign  of  their  slope.  Similar  results  have  been  observed 
for  other  values  of  the  maximum  allowable  acceleration  u“ax,  and  the  optimal  contour  can  be  expressed  more 


Fraction  of  face-sheet  mass  in  front  face  Fraction  of  face-sheet  mass  in  front  face 

%  /  (.m  +  >h)  *h  1  (7:  +  Vi) 

Fig.  17.  Contours  with  varying  mass  distribution  of  (a)  maximum  allowable  nondimensional  reflected  impulse  IR  and  (b)  maximum 
allowable  nondimensional  incident  impulse  /0  with  nondimensional  back-face  accelerations  limited  to  5™ax  =  5.  The  approximate 
treatment  of  FSI  is  not  applicable  in  the  hatched  region,  in  which  uFSi(f»)/iii(t»)  >  1.5  (PR=100,  P0=13,  pJp/,=  \'iQ, 
Paca/ \/  Poffp/8D  =  0.02206). 
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generally  as  a™1  =  a”ax-  This  means  that  for  a  given  i/0,  the  allowable  impulse  is  maximized  by  adjusting  the 
mass  distribution  so  that  the  nondimensional  acceleration  at  complete  densification,  a2nt,  equals  the  limiting 
value,  a”ax.  Noting  that  a2'1  =  1  /tj2,  this  optimal  mass  distribution  can  be  expressed  as  t]2  =  1  /a”ax,  which  sim¬ 
ply  corresponds  to  a  constant  mass  fraction  in  the  back  face.  For  a  given  ;/0,  this  criterion  specifies  a  mass 
fraction  in  the  front  face  of  rjl  =  1  —  r\Q  —  1  /a)"’.  which  decreases  with  increasing  i/0  and  reduces  to  zero  when 
t]Q  =  1  —  l/a”ax.  This  can  be  seen  in  Figs.  13a,  in  which  the  contour  for  =  5  intersects  the  vertical  axis  at 
t] o=l  —  1/5  =  0.8.  For  i/0  >1-1  / d"lax .  the  contours  of  Fig.  13a  indicate  that  the  allowable  impulse  is  max¬ 
imized  by  placing  all  of  the  face-sheet  mass  in  the  back  face  (i.e.,  setting  iq  —  0  and  r\2  =  1  —  r/0).  However,  it  is 
noted  that  as  i/0  becomes  large,  the  maximum  permissible  impulse  becomes  quite  insensitive  to  the  relative  dis¬ 
tribution  of  mass  between  the  front  and  back  faces.  For  ;/0  =  0.9,  for  example,  the  maximum  allowable 
impulse  for  rj\  /(iji  +  i]2)= 0.001  is  only  0.3%  larger  than  the  maximum  allowable  impulse  for  f/j/ 
0/1  +  m)  =  0.999. 


7.  Summary  and  conclusions 

This  paper  has  investigated  the  uniaxial  crushing  of  a  cellular  layer  sandwiched  between  solid  front  and 
back  faces,  with  air-blast  loading  applied  to  the  front  face  and  the  back-face  unrestrained.  While  previous 
studies  have  sought  to  minimize  deflections  for  sandwich  plates  with  face  sheets  of  equal  mass,  the  objective 
of  the  present  study  was  to  mitigate  shock  transmission  through  the  thickness  of  the  sandwich  plate,  and 
potential  advantages  of  unequal  face  sheets  were  explored  by  evaluating  response  characteristics  for  all  pos¬ 
sible  distributions  of  mass  between  the  front  and  back  faces  and  the  cellular  core.  To  investigate  the  crushing 
response  of  the  system,  an  analytical  model  was  developed  using  the  well-established  rigid,  perfectly-plastic, 
locking  (R-P-P-L)  idealization  for  the  cellular  core.  This  analytical  model  represents  the  propagation  of  a  pla¬ 
nar  densification  front  through  the  core  and  accounts  for  the  rate-dependent  energy  dissipation  associated 
with  the  densification  front,  which  can  be  many  times  larger  than  the  quasi-static  energy  dissipation.  Predic¬ 
tions  of  the  analytical  model  showed  excellent  agreement  with  computational  simulations  using  the  explicit 
finite  element  code  LS-DYNA,  while  a  fairly  small  value  of  the  material  viscosity  coefficient  was  required 
in  the  computations  to  avoid  additional  viscous  dissipation  due  to  spreading  of  the  densification  front. 

Fluid-structure  interaction  (FS1)  effects  were  treated  using  a  recent  result  due  to  Kambouchev  et  al.  (2006) 
for  air-blast  loading  on  solid  plates.  This  result  was  extended  to  sandwich  plates  by  defining  an  “equivalent” 
solid  plate  based  on  the  initial  acceleration  of  the  front  face.  This  FSI  approximation  neglects  increases  in  the 
effective  inertia  of  the  front  face  resulting  from  accretion  of  crushed  core  material  as  the  densification  front 
propagates,  and  it  was  found  to  be  inappropriate  when  the  mass  fraction  in  the  core  i/0  is  large  relative  to 
the  mass  in  the  front  face  i/i,  particularly  when  the  decay  period  of  the  pressure  pulse  is  comparable  to  the 
propagation  time  of  the  densification  front.  Within  its  domain  of  applicability,  this  approximation  relates 
the  nondimensional  reflected  impulse  7R  and  the  nondimensional  peak  reflected  pressure  PR,  which  are  applied 
to  the  front  face  of  the  sandwich  plate,  to  the  corresponding  values  I0  and  P0  associated  with  the  incident  air 
shock. 

The  capacity  of  the  sandwich  plate  to  mitigate  shock  transmission  is  limited  by  the  critical  impulse  required 
to  produce  complete  densification  of  the  core,  with  corresponding  nondimensional  reflected  and  incident  val¬ 
ues  denoted  7™‘  and  7™*.  For  7R  <  7™1,  the  cellular  core  serves  a  protective  function  by  limiting  the  stress  that 
reaches  the  back  face  to  the  plateau  stress  erP.  For  /R  >  /))“.  however,  reflection  of  the  densification  front  from 
the  back-face  results  in  large  stresses  that  can  exceed  the  peak  reflected  blast  pressure  by  several  times.  It  was 
observed  that  /R“  decreases  with  increasing  nondimensional  peak  pressure  PR,  being  smallest  in  the  limiting 
case  of  a  Dirac  delta  pulse  (PR  — >  oo),  which  corresponds  to  an  initial  velocity  imparted  to  the  front  face.  A 
limiting  nondimensional  peak  pressure  P1™  was  identified,  below  which  only  partial  densification  can  be 
achieved  even  if  the  pressure  is  maintained  indefinitely  f/j/1  — >  oo),  and  it  was  observed  that  densification 
of  the  core  cannot  initiate  if  PR  <  (1  —  ,  regardless  of  the  reflected  impulse.  Values  of  /Rnt  for 

PR  =  100  were  virtually  indistinguishable  from  those  for  the  Dirac  delta;  however,  for  large  ;/0,  the  rate-depen¬ 
dent  energy  dissipation  associated  with  densification  front  was  many  times  larger  for  the  Dirac  delta  pulse. 

The  influence  of  mass  distribution  on  the  critical  impulse  was  investigated,  and  it  was  found  that  both  7™1 
and  7™*  increase  with  increasing  mass  fraction  in  the  core  rj0.  For  a  given  ;/0,  both  and  7™'  were  found  to 
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increase  with  increasing  mass  fraction  in  the  front  face  i]  \ ,  with  corresponding  decreases  in  the  ratio 
^p°V^p°soiid  °f  work  done  by  the  pressure  pulse  on  the  sandwich  plate  to  the  work  done  on  a  correspond¬ 
ing  solid  plate  with  the  same  total  mass.  Impulse  reflection  coefficients  associated  with  the  critical  impulses 
were  observed  to  be  reduced  by  as  much  as  30%  relative  to  impulse  reflection  coefficients  for  a  rigid  wall, 
due  to  the  beneficial  effects  of  FSI.  In  spite  of  this,  the  ratio  W l°l /  ffip°‘olid  was  always  greater  than  unity,  indi¬ 
cating  that  impulse  reductions  due  to  FSI  are  insufficient  to  offset  the  enhanced  work  that  results  from 
increased  velocities  of  the  front  face.  The  ratio  Wl°l /Wl°\oM  is  maximized  for  small  values  of  rj0  and  ;/i,  with 
most  of  the  mass  in  the  back-face  (rj 2  — »  1),  for  which  values  of  ffip0t/^p°sohd  in  excess  of  10  are  observed,  and 
both  7™1  and  7™*  are  also  minimized  in  this  region.  Interestingly,  such  mass  distributions  correspond  precisely 
to  the  “blast  pendulum”  experiments  of  Hanssen  et  al.  (2002),  in  which  impulse  enhancement  was  observed.  It 
is  therefore  observed  that  such  configurations  are  particularly  well  suited  for  absorbing  energy  from  air  blasts 
and  particularly  poorly  suited  for  dissipating  energy  through  core  crushing,  which  helps  to  explain  their  poor 
performance  for  blast  mitigation. 

The  observed  increases  in  7™*  and  7™‘  with  corresponding  decreases  in  W{°1  /  W^oXlA  suggest  that  the  capac¬ 
ity  of  the  sandwich  plate  to  mitigate  shock  transmission  might  be  improved  by  increasing  the  mass  fraction  in 
the  front  face  i]\  and  reducing  the  mass  fraction  in  the  back-face  >]2-  However,  it  was  also  observed  that  reduc¬ 
ing  rj2  leads  to  increased  back-face  accelerations,  thus  decreasing  the  protection  provided  by  the  core  material. 
In  fact,  it  was  found  that  if  t]2  <  1/PR,  the  back-face  acceleration  can  exceed  the  peak  acceleration  of  the  center 
of  mass,  meaning  that  the  cellular  core  actually  amplifies  the  back-face  accelerations,  rather  than  reducing 
them.  Motivated  by  these  observations,  a  design  optimization  problem  was  posed  that  involves  finding  the 
mass  distribution  that  maximizes  the  nondimensional  impulse  that  can  be  absorbed  while  limiting  the  nondi- 
mensional  back-face  accelerations  to  a  specified  maximum  value,  denoted  a“ax.  In  the  domain  of  applicability 
of  the  FSI  approximation,  it  was  found  that  for  a  given  mass  fraction  in  the  core  i/0  (provided 
i'Iq  <  1  —  l/o“ax),  the  allowable  values  of  both  IR  and  70  are  maximized  by  setting  the  mass  fraction  in  the  back 
face  to  )]2  =  1  / a2  ,ix .  This  corresponds  to  a  mass  distribution  such  that  the  back-face  acceleration  at  complete 
densification  equals  the  maximum  allowable  value.  While  corresponding  values  of  70  could  not  be  evaluated 
for  large  ij0  due  to  the  inapplicability  of  the  FSI  approximation,  allowable  values  of  7R  were  found  to  be  quite 
insensitive  to  the  relative  distribution  of  mass  between  the  front  and  back  faces  for  rj0  >  1  —  l/a“ax.  More 
detailed  computational  modeling  of  FSI  effects  for  sandwich  plates  would  be  required  to  assess  whether  fur¬ 
ther  increases  in  the  allowable  values  of  I0  could  be  achieved  when  the  mass  fraction  in  the  front  face  /q  is 
small  relative  to  rj0. 
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